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ABSTRACT 

We simultaneously constrain cosmology and galaxy bias using measurements of galaxy 
■ abundances, galaxy clustering and galaxy-galaxy lensing taken from the Sloan Digi- 

tal Sky Survey. We use the conditional luminosity function (which describes the halo 
occupation statistics as function of galaxy luminosity) combined with the halo model 
(which describes the non-linear matter held in terms of its halo building blocks) to 
describe the galaxy-dark matter connection. We explicitly account for residual redshift 
space distortions in the projected galaxy-galaxy correlation functions, and marginalize 
over uncertainties in the scale dependence of the halo bias and the detailed structure of 
dark matter haloes. Under the assumption of a spatially flat, vanilla ACDM cosmology, 
we focus on constraining the matter density, fi m , and the normalization of the mat- 
ter power spectrum, erg, and we adopt WMAP7 priors for the spectral index, n s , the 
Hubble parameter, h, and the baryon density, fib- We obtain that fi m = 0.278^o;o26 
and as = 0.763tQ;o4g (95% CL). These results are robust to uncertainties in the ra- 
dial number density distribution of satellite galaxies, while allowing for non-Poisson 
satellite occupation distributions results in a slightly lower value for a 8 (0.744±g;g|f ). 
These constraints are in excellent agreement (at the la level) with the cosmic mi- 
crowave background constraints from WMAP. This demonstrates that the use of a 
realistic and accurate model for galaxy bias, down to the smallest non-linear scales 
currently observed in galaxy surveys, leads to results perfectly consistent with the 
vanilla ACDM cosmology. 

Key words: galaxies: halos — large-scale structure of Universe — dark matter — 
cosmological parameters — gravitational lensing — methods: statistical 



1 INTRODUCTION 

The last two decades have seen the emergence of a con- 
cordance cosmological model which describes the formation 
and evolution of cosmic structure in a scenario known as 
ACDM. In these cosmological models, gravity is described 
by General Relativity, dark matter and dark energy are the 
major constituents of the Universe (with normal 'baryonic' 
matter only contributing ~ 4 percent), and density pertur- 
bations are seeded by quantum fluctuations in a scalar field, 
the inflaton, that dominated the energy density shortly after 

* E-mail: cacciato@phys.huji.ac.il 



the Big Bang. In its most basic ('vanilla') form, the ACDM 
model assumes a flat geometry, dark energy is modeled as 
Einstein's cosmological constant, neutrino mass is assumed 
to be negligible, and the initial power spectrum of density 
perturbations is assumed to be a single power-law. Such 
ACDM cosmologies are described by 5 parameters: the en- 
ergy densities (in terms of the critical density) of baryons, 
fib, and cold dark matter, fidm, the spectral index, n s , and 
normalization, as, of the initial power spectrum, and the 
Hubble parameter, h = H o /(100 kms" 1 Mpc" 1 ). The flat 
geometry implies that fiA = 1 — fi m = 1 — fib — fidm, and is 
strongly supported by the location of the first peak in the 
angular power spectrum of cosmic microwave background 
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(CMB) temperature fluctuations (e.g., Balbi et al. 2000; 
Lange et al. 2001; Pryke et al. 2002; Netterfield et al. 2002; 
Ruhl et al. 2003) combined with the results on the Hub- 
ble constant from the Hubble Key Project (Freedman et 
al. 2001). 

All these cosmological parameters have now been con- 
strained at the few (~ 3-20) percent level by a variety of 
probes, including, among others, temperature anisotropies 
in the CMB (e.g., Spergel et al. 2003, 2007; Reichardt et 
al. 2009; Dunkley et al. 2009; Komatsu et al. 2009, 2011), 
Cepheid distances (e.g., Freedman et al. 2001; Sandage et 
al. 2006; van Leeuwen et al. 2007), high redshift super- 
novae la (e.g., Riess et al. 1998; Pcrlmutter et al. 1999, 
Astier et al. 2006; Kowalski et al. 2008), measurements of 
the primordial deuterium abundance (e.g., Buries, Nollett & 
Turner 2001; O'Meara et al. 2006), cluster abundances (e.g., 
Vikhlinin et al. 2009; Rozo et al. 2010; Benson et al. 2011), 
cosmic shear (e.g., Benjamin et al. 2007; Fu et al. 2008; Lin 
et al. 2011), the integrated Sachs-Wolfe effect (e.g., Ho et 
al. 2008; Giannantonio et al. 2008), the Lya forest (e.g., Viel, 
Weller & Haehnelt 2004; McDonald et al. 2005; Desjacques 
& Nusser 2005), and strong gravitational lensing (e.g., Koop- 
mans et al. 2003; Oguri et al. 2008). Despite some tension 
between a few subsets of all these independent constraints 
(see e.g., Dunkley et al. 2009 for a comprehensive overview), 
overall they are in good mutual agreement, giving rise to the 
notion of a true concordance cosmology. 

Another potentially powerful probe for cosmology is the 
(large-scale) distribution of galaxies. Although stars make 
a negligible contribution to the total energy density of the 
Universe, the light from stars in galaxies can be observed di- 
rectly and over cosmological scales, making galaxies useful 
tracers of the underlying dark matter density field. Unfortu- 
nately, the connection between galaxies and (dark) matter 
is muddled by the fact that galaxies are biased tracers of the 
mass distribution. The main problem is that this 'galaxy 
bias' is known to be extremely complicated: it is stochastic 
(e.g., Dekel & Lahav 1999; Tegmark & Bromley 1999), de- 
pends on galaxy properties such as luminosity, color and/or 
morphological type (e.g., Park et al. 1994; Guzzo et al. 2000; 
Norberg et al. 2001, 2002; Zehavi et al. 2005, 2011; Wang 
et al. 2007), and is scale dependent on small scales (e.g., 
Pervical et al. 2007; Reid, Spergel & Bode 2009; Cacciato 
et al. 2012). Based on these considerations, it is not surpris- 
ing that galaxy bias is generally considered a nuisance when 
using galaxies to constrain cosmology. However, galaxy bias 
also contains a wealth of information regarding galaxy for- 
mation, especially on small scales (e.g., Cacciato et al. 2012). 
After all, it is the (poorly understood) physics of galaxy for- 
mation that determines where, how and with what efficiency 
galaxies form within the dark matter density field (see Mo, 
van den Bosch & White 2010). Therefore, ideally, one would 
like to simultaneously solve for cosmology and galaxy bias. 

Early attempts to do so used galaxy power spectra mea- 
sured from the two-Degree Field Galaxy Redshift Survey 
(2dFGRS; Colless et al. 2003) and/or Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000) combined with a phenomeno- 
logical fitting function for the non-linear, scale-dependent 
bias and marginalized over its free parameters (e.g., Cole et 
al. 2005; Tegmark et al. 2006; Padmanabhan et al. 2007). 
However, this approach has two problems. First, the fitting 
function used to describe the galaxy bias has no logical con- 



nection to galaxy formation. For this reason the free pa- 
rameters used to parameterize the galaxy bias are merely 
considered as nuisance parameters. Second, as pointed out 
by Sanchez & Cole (2008) and Yoo et al. (2009), the par- 
ticular fitting function used is often a poor description of 
the true scale dependence of galaxy bias, causing biased es- 
timates of the cosmological parameters, especially for non- 
vanilla ACDM cosmologies (e.g., Hamann et al. 2008). It 
is generally believed that the latter problem explains why 
there is still some tension between cosmological parameters 
(mainly f2 m ) inferred from different galaxy power spectra, 
or from using the same power spectrum, but using data cov- 
ering different scales (e.g., Percival et al. 2007; Sanchez & 
Cole 2008; Dunkley et al. 2009) 

This indicates that it is prudent to use a realistic, physi- 
cally motivated model for galaxy bias. Under the assumption 
that all galaxies reside in dark matter haloes, a natural, re- 
alistic model for galaxy bias is provided by halo occupation 
models, which describe, in a statistical sense, how galax- 
ies are distributed over dark matter haloes (e.g., Jing, Mo & 
Borner 1998; Peacock & Smith 2000; Scoccimarro et al. 2001; 
Berlind & Weinberg 2002; Yang, Mo & van den Bosch 2003). 
When combined with the halo model, which describes the 
(non-linear) matter distribution in terms of the dark mat- 
ter halo building blocks (e.g., Neyman & Scott 1952; Seljak 
2000; Ma & Fry 2000; Cooray & Sheth 2002), these halo 
occupation models provide a complete, accurate, and easy- 
to-interpret description of galaxy bias, all the way from the 
large, linear scales down to the small, non-linear scales of 
individual dark matter haloes. 

The problem with this approach, though, is that the 
halo occupation models are also cosmology dependent, so 
that one typically needs constraints in addition to just large 
scale clustering. Several attempts have been made along 
these lines. Abazajian et al. (2005) simply used strong priors 
from the WMAP results. Yang et al. (2004) and Tinker et 
al. (2007) used peculiar velocities as inferred from the red- 
shift space distortions in the two-point correlation function, 
and argued for a relatively low value of as ~ 0.7 — 0.75 (for 
a vanilla ACDM cosmology with f2 m = 0.3). Similar con- 
clusions were reached by van den Bosch et al. (2003b) and 
Tinker et al. (2005), who used constraints on the (average) 
mass-to-light ratios of clusters, rather than peculiar veloci- 
ties. Somewhat puzzling, a more recent analysis by Tinker et 
al. (2012) used the mass-to-number ratio of clusters, rather 
than the mass-to-light ratios, and found a relatively high 
value for as of 0.85 ± 0.05 (again for f2 m = 0.3). 

In this paper we use a combination of galaxy clustering 
and galaxy-galaxy lensing, as well as constraints on galaxy 
abundances, in order to simultaneously constrain cosmology 
and galaxy bias. Since galaxy-galaxy lensing probes the mass 
associated with the lensing galaxies, this is similar to us- 
ing mass-to-light ratios as constraints. It has the advantage, 
though, that it probes mass-to-light ratios over a wide range 
in halo masses, and that the same halo occupation model 
used to compute the clustering of galaxies can also be used 
to compute the galaxy-galaxy lensing signal (e.g., Guzik & 
Seljak 2002; Yoo et al. 2006; Cacciato et al. 2009). A first ap- 
plication of this idea by Seljak et al. (2005) , when combined 
with WMAP constraints, yielded a relative high value for 
a s of 0.88 ± 0.06 (for fi m ~ 0.25). More recently, two differ- 
ent analyses based on the same galaxy-galaxy lensing data 
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by Cacciato et al. (2009) and Li et al. (2009) both argued 
that a flat ACDM cosmology with (fi m , a 8 ) = (0.238, 0.734) 
is in much better agreement with the data than a (0.3, 0.9) 
model, thus favoring again a relatively low value for as- In 
this paper we improve on all these previous methods by (i) 
simultaneously constraining galaxy bias and cosmology, (ii) 
using a much more accurate analytical model, (iii) using the 
latest clustering data, (iv) modelling the latter accounting 
for residual redshift space distortions, (v) marginalizing over 
uncertainties related to the detailed structure of dark mat- 
ter haloes and the scale dependence of the halo bias, and 
(vi) subjecting the analysis to a number of detailed tests 
that address how the results depend on certain assumptions 
inherent to the model. 

This paper is the third in a series. In van den Bosch 
et al. (2012; hereafter Paper I), we presented the analytical 
model, which we calibrated and tested using detailed mock 
catalogs constructed from high resolution JV-body simula- 
tions. We demonstrated that our analytical model is accu- 
rate at the level of 10% or better over the entire range of 
scales covered by the data. In More et al. (2012b; hereafter 
Paper II), we presented a Fisher matrix analysis to iden- 
tify parameter-degeneracies and to assess the accuracy with 
which various cosmological parameters can be constrained 
using our methodology. We demonstrated that the method 
can simultaneously constrain halo occupation statistics and 
cosmology, and we forecasted that, using existing data from 
the SDSS, we should be able to put constraints on f2 m and 
as that are among the tightest ever achieved. In this pa- 
per we apply our method to existing data from the SDSS^. 
Although, as demonstrated in paper II, our method is also 
able to constrain extensions to the vanilla ACDM cosmol- 
ogy, such as neutrino mass and the equation of state of dark 
energy, in this paper we focus solely on vanilla ACDM cos- 
mologies, and in particular on constraining the combination 
(Qm_,as). We defer constraining neutrino mass, dark energy, 
and other modifications to the vanilla ACDM cosmology to 
future papers. 

This paper is organized as follows. In §2, we introduce 
the SDSS data used to constrain our models. In §3, we briefly 
review our analytical model to compute the galaxy luminos- 
ity function, the galaxy-galaxy correlation function, and the 
galaxy-galaxy lensing signal using the halo model combined 
with a model that describes halo occupation statistics as 
function of galaxy luminosity. The Bayesian analysis, used to 
infer posterior distributions for the cosmological parameters 
and for the parameters that describe the halo occupation 
statistics, is described in §4. Our main results are presented 
in §5, while §6 describes a number of tests that address the 
sensitivity of our results to several model assumptions. We 
summarize our findings in §7. 

Throughout this paper, unless specifically stated other- 
wise, all radii and densities are in comoving units, and log 
is used to refer to the 10-based logarithm. Quantities that 
depend on the Hubble parameter are written in units of h, 
defined above. 



2 DATA 

The data used to constrain our models consists of three com- 
ponents: galaxy abundances, in the form of the galaxy lumi- 
nosity function, galaxy clustering, in the form of projected 
correlation functions for six different luminosity bins, and 
galaxy-galaxy lensing, in the form of excess surface densi- 
ties (ESD), once again for six different luminosity bins. All 
these measurements are obtained from the Sloan Digital Sky 
Survey (SDSS; York et al. 2000). 

We use the r-band galaxy luminosity function, 5>(L), 
(hereafter LF) of Blanton et al. (2003a)*, sampled at 32 
magnitudes in the range —23.0 < 0,1 M r — 5 log h < —18.0, 
where 0,1 M r indicates the r-band magnitude of galaxies 
K+E corrected to z = 0.1 following the procedure of Blanton 
et al. (2003b). For each of these magnitude bins we use the 
(statistical) errors on $ quoted by Blanton et al. (2003a). 
Unfortunately, we do not have a full covariance matrix for 
this data set. There are two main sources of covariance 
for the luminosity function: magnitude errors, which cause 
covariance between neighboring magnitude bins, and large 
scale structure ('sample variance')^. The LF of Blanton et 
al. is sampled at magnitude intervals of 8m r = 0.00986 mag, 
which is much smaller than the typical error on individual 
magnitudes (~ 0.05 mag). In order to suppress the covari- 
ance due to these magnitude errors, we only sample the LF 
of Blanton et al. at 32 magnitude intervals of 8m r = 0.1578, 
roughly three times as large as the typical magnitude error. 
This should eliminate virtually all covariance due to errors 
in the magnitudes of individual galaxies. This leaves the co- 
variance due to sample variance, which can be 'modeled' 
to reasonable accuracy as an up or down shift of the en- 
tire LF (i.e.,. the error bars are fully covariant, e.g. Blanton 
et al. 2003a). Throughout this study, we ignore this sample 
variance, which effectively implies that we assume that the 
volume probed by the SDSS is a fair representation of the 
average Universe (see discussion in §6.4). 

The galaxy clustering data used in this paper are taken 
from Zehavi et al. (2011), and based on the SDSS DR7 
(Abazajian et al. 2009). Using all galaxies in the main galaxy 
sample with apparent magnitudes 14.5 < m r < 17.6, Zehavi 
et al. measured the projected correlation functions, w p (r p ), 
over the radial range 0.2/i _1 Mpc < r p < 20h~ 1 Mpc for 
six volume limited samples (see Table 1). These have been 
obtained according to 

/''"max 

w p (r p ) = 2 Cg g (r p ,r 7r )dr 7r . (1) 

Jo 

with r max = 40/i -1 Mpc or 60/i _1 Mpc, depending on the lu- 
minosity sample used (see Table 1). Here r p is the projected 
separation between two galaxies, tv is the redshift-spacc 
separation along the line-of-sight, and £ gg (r p , tv) is the two- 
dimensional correlation function, which is anisotropic due to 
the presence of peculiar velocities. As discussed in Paper I, 
the fact that r max is finite results in residual redshift space 
distortions that need to be corrected for in the modeling 
(see also § 6.3). The errors on w p (r p ) are characterized by 



t A preliminary version of the main results presented in this 
paper were published in a conference proceedings by More et 
al. (2012a). 



t available at http://cosmo.nyu.edu/mbl44/lf.html 

§ Following Scott, Srednicki & White (1994) we use the term 

'sample variance', rather than the more common 'cosmic variance' 
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Table 1. Properties of the projected correlation function data 



1 M r - 5 log h 




% max 


(z) 




^raax 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(-18,-17] 


0.011 


0.026 


0.021 


5,965 


40 


(-19,-18] 
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18,200 


40 


(-20,-19] 


0.027 


0.064 


0.050 


44,348 


60 


(—21, -20] 


0.042 


0.106 


0.082 


108,629 


60 


(-22,-21] 


0.066 


0.159 


0.123 


73,746 


60 


(-23, -22] 


0.103 


0.245 


0.187 


10,251 


60 



The volume-limited galaxy samples used by Zehavi et al. (2011) 
to measure the projected correlation function, w p (r p ). For each of 
these samples column (1) lists the magnitude range in 0,1 M r — 
51og/i, columns (2), (3) and (4) list the minimum, maximum 
and mean rcdshifts, respectively, column (5) lists the number of 
galaxies, and column (6) lists the value of r max (in ft -1 Mpc), 
defined in Eq. (1). 



their full covariance matrices, kindly provided to us in elec- 
tronic format by I. Zehavi, and obtained from 144 spatially 
contiguous subsamples using the jackknife technique. 

Finally, for the galaxy-galaxy lensing data we use the 
excess surface densities, AT,(R), covering the radial range 
0.04ft _1 Mpc < R < 2\r x Mpc, obtained by Seljak et 
al. (2005) and Mandelbaum et al. (2006), and kindly pro- 
vided to us in electronic format by R. Mandelbaum. These 
measurements have been obtained using a catalogue of 
351, 507 lens galaxies with apparent magnitude m r < 17.77 
in the redshift range 0.02 < z < 0.35 taken from the main 
galaxy catalogue of the SDSS Data Release 4 (Adelman- 
McCarthy et al. 2006). This sample is split in 8 (flux- limited) 
luminosity bins. However, since both the faintest and bright- 
est bins have extremely poor signal-to-noise, we only use 
the 6 intermediate luminosity bins listed in Table 2. We re- 
fer the reader to Mandelbaum et al. (2006) for a detailed 
description of the data and of the method used to deter- 
mine the excess surface density profiles. Since the covariance 
in AE(J?) is only very small over the radial scales covered 
by the data (R. Mandelbaum, private communication), we 
only use the diagonal elements in what follows. Throughout 
this paper we refer to these three data sets as the LF, WP 
and ESD data. Note that each of these data sets adopted 
a flat ACDM cosmology with Q m = 0.3 (LF and WP) or 
f2 m = 0.27 (ESD) when computing distances and/or abso- 
lute magnitudes. Changing the assumed cosmology, as we 
do in our analysis, in principle therefore also has a small 
impact on the observational measurements by changing the 
distance-redshift relation and thus shifting galaxies among 
luminosity bins and galaxy pairs among radial separation 
bins. However, in this paper we restrict ourselves to cos- 
mologies that only differ mildly from a ACDM cosmology 
with fi m = 0.3. Even at our outer redshift limit of z = 0.2, 
the effect of lowering fi m from 0.3 to 0.25 is only 1 percent 
in distance, so the measurements used here are effectively 
independent of cosmological parameters within their obser- 
vational uncertainties. Therefore, we ignore this small effect 
in what follows, and always use the observational data as 
described above. 



Table 2. Properties of the excess surface density data 



0-1 A/f 


5 log h 


1 -A 


AT 

1 * lens 


(i) 




(2) 


(3) 


(-19.0,- 


-18.0] 


0.047 


29,730 


(-20.0, - 


-19.0] 


0.071 


85,766 


(-21.0, - 


-20.0] 


0.10 


141,976 


(-21.5, - 


-21.0] 


0.14 


60,994 


(-22.0, - 


-21.5] 


0.17 


34,920 


(-22.5, - 


-22.0] 
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The flux-limited galaxy samples used by Mandelbaum ct 
al. (2006) to measure the excess surface density profiles, AE(R). 
For each of these samples column (1) lists the magnitude range in 
01 M r — 51og/i, column (2) lists the mean redshift, and column 
(3) lists the number of lens galaxies used. 

3 MODEL DESCRIPTION 

As described in detail in Paper I, the observables^ &(L,z), 
w p (r p ,z) and AS(_R, z) can be computed for a given cos- 
mology, which determines the properties of the dark mat- 
ter distribution (e.g, halo mass function, halo bias function, 
halo density profiles), and a given description of the galaxy- 
dark matter connection (i.e., halo occupation statistics). The 
dependence on redshift is included to emphasize that we 
model each observable at its mean redshift. For the lumi- 
nosity function we adopt z = 0.1, while the mean redshifts 
for the different luminosity bins of the w p (r p ) and AE(i?) 
measurements are listed in Tables 1 and 2, respectively. 

In this section we give a concise overview of our model. 
Readers interested in a more thorough description are re- 
ferred to Paper I, while those readers that are already fa- 
miliar with our model, or that are mainly interested in the 
results, may want to skip this section and proceed immedi- 
ately to § 4. 

3.1 Cosmological Parameters 

Throughout this paper we consider 'vanilla' ACDM cosmolo- 
gies in which gravity is described by standard General Rel- 
ativity, neutrino mass is negligible, the initial power spec- 
trum is a single power-law, and dark energy is modeled as 
Einstein's cosmological constant with J1a = 1 — fi m - These 
cosmologies are completely specified by five parameters; the 
z — matter density in units of the critical density, f2 m , 
the normalization of the matter power spectrum, as, the 
Hubble parameter h — _ffo/(100 kms -1 Mpc -1 ), the initial 
spectral index of the matter power spectrum, n B , and the 
z = baryon density in units of the critical density, fib- 
Hence, our cosmological model parameters are described by 
the vector 

A coenM = (n m ,ff 8 ,fc,n.,n b fc 2 ). (2) 

Note that the baryon density only enters in our analysis in 
the transfer function, and always in the combination fib h 2 . 

^ The projected correlation function and the excess surface den- 
sity are modeled per luminosity bin [Li, Li]- However, to keep the 
notation concise, we do not explicitly write down the dependence 
on L\ and Li. 
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The main goal of this paper is to use the observational 
data on w p (r p ) and AE(7?) discussed in §2 to con- 

strain O m and as- Throughout this paper we will there- 
fore not use any priors on these two parameters. For h, n B 
and f2b, on the other hand, we include prior information 
from the seven year analysis of the cosmic microwave back- 
ground data from WMAP (hereafter WMAP7; Komatsu et 
al. 2011), as described in §4 below. For this reason we refer 
to h, n B and fi b as our secondary cosmological parameters 
in what follows. 



3.2 The Conditional Luminosity Function 

Under the assumption that each galaxy resides in a dark 
matter halo, $(L,z), w p (r p ,z) and AE(P, z) can be com- 
puted using a statistical description of how galaxies are dis- 
tributed over dark matter haloes of different mass. To that 
extent we use the conditional luminosity function (here- 
after CLF) introduced by Yang et al. (2003). The CLF, 
<&(L\M)dL, specifies the average number of galaxies with 
luminosities in the range L ± dL/2 that reside in a halo of 
mass M. We split the CLF in two components, 



$(L|M) = $ C (L|M) + $ S (L|M) 



(3) 



where 3> C (L|M) describes the contribution due to central 
galaxies (defined as those galaxies that reside at the cen- 
ter of their host halo), while & S (L\M) characterizes satellite 
galaxies (those that orbit around a central). Throughout we 
ignore a potential redshift dependence of the CLF. Since the 
data that we use to constrain the CLF only covers a nar- 
row range in redshift, this assumption will not have a strong 
impact on our results. 

Our parameterization of the CLF model is motivated by 
the results obtained by Yang, Mo & van den Bosch (2008) 
from a large galaxy group catalogue (Yang et al. 2007) 
extracted from the SDSS Data Release 4, and by Tal et 
al. (2012) from a study of the luminosity function of satel- 
lite galaxies of luminous red galaxies. In particular, the CLF 
of central galaxies is modeled as a log-normal function: 



$c(L\M) dL= 1 °i 6 exp 



2-7T <T C 



(logL-logLc) 2 



2ai 



dL 

L ' 



(4) 



and the satellite term as a modified Schechter function: 



$ s (L|M)dL = 



exp 



which decreases faster than a Schechter function at the 
bright end. Note that L c , a c , cf>*, a s and L* are all func- 
tions of the halo mass M. 

Following Cacciato et al. (2009), and motivated by the 
results of Yang et al. (2008) and More et al. (2009, 2011), 
we assume that a c , which expresses the scatter in logL of 
central galaxies at fixed halo mass, is a constant (i.e. is in- 
dependent of halo mass and redshift). In addition, for L c , 
we adopt the following parameterization; 



L C (M) 



(M/Mi) 71 



(6) 



[l-KM/Mi)] 71 - 72 

Hence, L c oc M 71 for M « Mi and L c oc M 72 for 
M 3> M\. Here Mi is a characteristic mass scale, and 
Lo = 2 71 ~ 72 L c (Mi) is a normalization. 
For the satellite galaxies we adopt 



L*(M) = 0.562L C (M) . 
q s (M) = a s 



(7) 
(8) 



(i.e., the faint-end slope of <E> S (L|M) is independent of mass 
and redshift), and 



log[0:(M)] = b + &i (log Mia) + MlogMi 2 ) 2 , 



(9) 



with M12 = M/(10 12 /i _1 M Q ). Note that neither of these 
functional forms has a physical motivation; they merely were 
found to adequately describe the results obtained by Yang 
et al. (2008) from the SDSS galaxy group catalog. 

To summarize, our parameterization of the CLF thus 
has a total of nine free parameters, characterized by the 
vector 



A 



(log Mi, log Lo,7i,72,<7c,as,&o,&i,&2) • 



(10) 



3.3 Galaxy Luminosity Function 



Once the CLF is specified, the galaxy luminosity function 
at redshift z, $(L, z), simply follows from integrating over 
the halo mass function, n(M, z); 



I 



&{L,z)= / $(L\M)n(M,z)dM . 



(11) 



In what follows, we will always be concerned with galaxies in 
a specific luminosity interval [Li , L2] . The average number 
density of such galaxies follows from the CLF according to 



fig(z) = J (N s \M)n(M,z)dM , 
where 

(N S \M) = I $(L|M)dL, 



J L-i 



(12) 



(13) 



is the average number of galaxies with Li < L < L2 that 
reside in a halo of mass M. 



3.4 Power-Spectra 

The first step towards computing the projected correlation 
functions, w p (r p , z), and ESD profiles, AE(7?, z), is to com- 
pute the galaxy-galaxy power spectrum, P gg (fc, 2), and the 
galaxy-matter cross power spectrum, P gm (fc, z). These power 
spectra are the Fourier space analogs of the galaxy-galaxy 
and galaxy-matter correlation functions. 

The galaxy-galaxy power spectrum can be expressed as 
a sum of the one-halo (lh) and the two-halo (2h) terms, each 
of which can be further subdivided based upon the type of 
galaxies (central or satellite) that contribute to the power 
spectrum, i.e., 

P gg (fc) = 2P c 1 s h (fc) + P s 1 s h (fc) 

+ P < ? c h (fc) + 2P < 2 3 h (fc) + P s 2 s h (fc). 



(14) 



Similarly, the galaxy-matter power spectrum can be written 

as 

P gm (fc) = Plt(k) + P£(k) + Pci(fe) + PS(fc) . (15) 

As shown in paper I, these terms can be written in compact 
form as 
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P xy (k, z 



I 



,z) = J n, 

-I 



(k,M,z)H y (k,M,z)n(M,z)AM, (16) 



dM 1 H x (k,M 1 ,z)n(M 1 ,z) 



AM 2 H y (k, M 2 , z) n(M 2 ,z) Q(k\M 1 ,M 2 , z) , (17) 



where 'x' and 'y' are either 'c' (for central), 's' (for satellite), 
or 'm' (for matter), Q(k\Mi,M 2 ,z) describes the power- 
spectrum of haloes of mass M\ and M 2 (see Appendix A), 
and we have defined 



M 

Hm(k, M,z) = — w h (fc|M, z) 

pm 



Hc(k,M,z)=H c (M,z) 



(N C \M) 
n e (z) 



(18) 
(19) 

(20) 



and 

H s (k, M, z) = {Ns \ M) u s (k\M, z) . 
n g (z) 

Here (N C \M) and (N S \M) are the average number of central 
and satellite galaxies in a halo of mass M, which follow from 
Eq. (13) upon replacing $(L|M) by $ c (LjM) and $ S (L|M), 
respectively. Furthermore, u s (k\M) is the Fourier transform 
of the normalized number density distribution of satellite 
galaxies that reside in a halo of mass M, and u h (k\M) is 
the Fourier transform of the normalized density distribution 
of dark matter within a halo of mass M. 



3.5 Computing w p (r p ,z) and AE(ii, z) 

Once P ss (k,z) and P sul (k,z) have been determined, it is 
fairly straightforward to compute the projected galaxy- 
galaxy correlation function, w p (r p ,z), and the excess sur- 
face density (ESD) profile, AT,(R, z). We start by Fourier 
transforming the power-spectra to obtain the two-point cor- 
relation functions: 



(21) 



where 'x' is either 'g' (for galaxies) or 'm' (for matter). 

The projected galaxy-galaxy correlation function, 
w p (R,z), is related to the real-space galaxy-galaxy corre- 
lation function, £g g (r, z), according to 



w P (r p ,z) = 2 / £ gg (r p ,r 7r ,z)Ar n 



= 2 



2 



&i(s,z)V 2 i(r 7r /s) drv. 



(22) 



Here r max is the maximum integration range used for the 
data (see Table 1), rv is the separation perpendicular to the 
line-of-sight, s = r p + r% is the separation between the 
galaxies in redshift space, Vi(x) is the Z th Legendre polyno- 
mial, and £o, £2, and £4 are given by 



Zo(r,z)= (l + |/3+ |/3 2 ) Znfaz), (23) 
6(r, z) = (J/3 + K„(r, z) - 3J 3 (r, z)] , (24) 



where 



1 ^ ^ 

^ g (r,z) + ^.Mr,z)-^Mr,z) , (25) 



J "( r ' 2 ) = ^"/ ^{V,z)y n Ay. 



and 



/3 = /3(z) = 



1 /dlnD\ ^ fim 6 (z) 
Uz) \ dlna ) z ~ b(z) 



(26) 



(27) 



with a = 1/(1 + z) the scale factor, D(z) the linear growth 
rate, and 



Hz) 



n g (z) 



(N g \M) b h (M, z) n(M, z) AM , 



(28) 



the mean bias of the galaxies in consideration. Note that 
Eq. (22) accounts for the large scale redshift space distor- 
tions due to infall (the 'Kaiser'-effect), which is necessary 
because the measurements for w p (p) have been obtained us- 
ing a finite r max . Note that the £gg( r ) in Eqs. (23)-(25) is the 
non-linear galaxy-galaxy correlation function. Although the 
Kaiser formalism (Kaiser 1987) is only strictly valid in the 
linear regime, this simple modification results in a more ac- 
curate correction for the residual redshift space distortions 
(see Paper I for details). As shown in Paper I, not taking 
these residual redshift space distortions into account results 
in systematic errors that can easily exceed 20 percent on 
large scales (r p > lO/i -1 Mpc), causing systematic errors in 
the inferred galaxy bias (see also Padmanabhan, White & 
Eisenstein 2007; Norberg et al. 2009; Baldauf et al. 2010; 
More 2011). 

Finally, the excess surface density profile, AH(R,z), is 
defined as 



f S(i?', z) R' AR' — E(-R, z) 
Jo 



(29) 



Here S(i?, z) is the projected surface mass density, which 
is related to the galaxy-dark matter cross correlation, 
£ gm (r, z), according to 



X(R,z) = 2p r 



■f 

Jr. 



[l + £ gm (r,z)] 



rdr 



Vr 2 - R 2 



(30) 



3.6 Model Ingredients 

In this subsection we briefly describe the ingredients of our 
model. A more detailed description can be found in Paper I. 

In our fiducial model we require both the linear and 
the non-linear power spectra of the matter distribution, 
PuSn{k,z) and P mm (fc, 2), respectively. The latter enters in 
the computation of the halo-halo correlation function, as 
described in detail in Appendix A. Throughout we compute 

Pmm(fc, z) using the fitting formula of Smith et al. (2003) I , 
while for i*(it,«) we use the linear transfer function of 
Eisenstein & Hu (1998), neglecting any contribution from 
neutrinos and assuming a CMB temperature of 2.725K 
(Mather et al. 1999). 



II We use the small modification suggested on John Peacock's 
website http://www.roe.ac.uk/~jap/halocs/, although we have 
verified that this has no significant impact on any of our results. 
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Table 3. List of models 



Model 


v> 


V 


■R 


V 


RRSD 


SGW 






Xesd 


X 2 ■ 

^prior 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


Fiducial 


5(0.9,0.15) 


5(0.0,0.1) 


1.0 


1.0 


yes 


yes 


47 


236 


117 


4 


2.36 


Radial 


£(0.9,0.15) 


5(0.0,0.1) 


5(1.0,0.25) 


1.0 


yes 


yes 


40 


240 


118 


5 


2.36 


Poisson 


5(0.9,0.15) 


5(0.0,0.1) 


1.0 


5(1.0,0.1) 


yes 


yes 


52 


229 


113 


7 


2.35 


noRSSD 


5(0.9,0.15) 


5(0.0,0.1) 


1.0 


1.0 


no 


yes 


48 


218 


114 


5 


2.25 


noSGW 


5(0.9,0.15) 


5(0.0,0.1) 


1.0 


1.0 


yes 


no 


32 


155 


111 


16 


1.99 



List of models described in the text. Column (1) indicates the model name. Columns (2) to (5) indicate the priors used for the parameters 
i/>, n, TZ and Ap. Here Q(x,a) indicates a Gaussian prior centered on x and with a standard deviation a. Column (6) indicates whether 
or not the model corrected for residual redshift space distortions (RRSD), and Column (7) indicates whether or not the Sloan great wall 
(SGW) is included in the clustering data. Finally, Columns (8) to (12) list the various x 2 values of the best-fit models in each MCMC 
chain. The reduced x 2 value in Column (12) is computed for N^ { = 171 degrees of freedom (see §5). 



Throughout, we define dark matter haloes as spheres 
with an average density of 200 times the background density. 
We assume that their density profiles follow the NFW profile 
(Navarro, Frcnk & White 1997), with a concentration- mass 
relation given by 



c(M,z) = (1 + 7?) c{M,z). 



(31) 



Here c(M, z) is the average concentration-mass relation of 
Maccio et al. (2007), properly converted to our definition of 
halo mass. We treat n as a free nuisance parameter which 
accounts for (i) the fact that there is an uncertainty of ~ 10 
percent in the average concentration mass relation as ob- 
tained by different authors (e.g., Navarro, Frenk & White 
1997; Eke, Navarro & Steinmetz 2001; Bullock et al. 2001; 
Maccio et al. 2007; Zhao et al. 2009), and (ii) realistic dark 
matter haloes are triaxial, rather than spherical, have sub- 
structure, and have scatter in the concentration-mass rela- 
tion. As shown in Paper I, setting rj = ±0.2 modifies the 
one-halo term of £ gm (r) by more than 20 percent on small 
scales (r < O.lfe -1 Mpc). As discussed at length in Paper I, 
an uncertainty in ij of this amount is more than adequate 
to capture the inaccuracies in our model that arise from the 
various oversimplifications and uncertainties regarding the 
structure of dark matter haloes. Hence, in what follows, we 
adopt a Gaussian prior on rj, centered on rj — and with a 
standard deviation a v — 0.1. 

For the halo mass function, n(M,z), and the halo bias 
function, bh(M, z), we use the fitting functions of Tinker et 
al. (2010). Note that these functions obey the normalization 
condition 



Mn(M,z)b h (M,z) AM = 1 



(32) 



as required by the fact that, on average, on large scales 
matter is unbiased with respect to itself. As described in 
Appendix A, the scale dependence of the halo bias in the 
quasi-linear regime is described by a modified version of the 
empirical fitting function of Tinker et al. (2005) . This modi- 
fication is needed to account for the fact that we use a differ- 
ent definition of dark matter haloes, and is characterized by 
one free 'nuisance' parameter, ip. In Paper I we calibrated 
ip using numerical iV-body simulations of structure forma- 
tion in a ACDM cosmology, and found that ip = 0.9 can 
accurately fit the simulation results. In order to account for 
the fact that we cannot rule out that ip is cosmology de- 



pendent, we include an uncertainty of ~ 15 percent on ip in 
our cosmological analysis. In particular, we treat ip as a free 
parameter, but adopt a Gaussian prior centered on ip — 0.9 
and with = 0.15. 

For our fiducial model, we assume that the radial num- 
ber density distribution of satellite galaxies follows that of 
the dark matter particles, i.e., u s (r\M) = Uh(r\M). In other 
words, we assume that satellite galaxies follow a NFW pro- 
file with the same concentration- mass relation (Eq. [31]) as 
dark matter haloes. In §6.1 we relax this assumption and 
examine how changes in u s (r\M) impact on our results. 

Finally, we emphasize that the expression for the 1-halo 
term of the galaxy-galaxy correlation function (Eq. [16]) has 
made the implicit assumption that the halo occupation num- 
ber of satellite galaxies obeys Poisson statistics, i.e., that 
(N S (N S - \)\M) = (N S \M) 2 (see Paper I for details). In §6.2 
we will relax this assumption, and explore how deviations of 
P(N B \M) from a Poisson distribution impact on our results. 



4 METHODOLOGY 

The main goal of this paper is to obtain constraints on the 
cosmological parameters f2 m and as, and the halo occupa- 
tion distribution, as characterized by the CLF, using the 
SDSS data described in §2 and the analytical model de- 
scribed in §3. We use Bayesian inference techniques to de- 
termine the posterior probability distribution P(\\T>) of the 
model parameters X, given the data V. According to Bayes' 
theorem, 



P(X\V) 



P{V) 



(33) 



where P(V\X) is the likelihood of the data given the model 
parameters, P(X) is the prior probability of these parameters, 
and 



/ 



P(p)= / P(V\X)P(X)dX, 



(34) 



is the marginal probability of the data, also called evidence 
for the model. Since, we do not intend to perform model 
selection, the evidence just acts as a normalization constant 
which need not be calculated. Therefore the posterior dis- 
tribution P(X\T>) is given by 
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-18 



-20 



-22 



o.i 



M r - 5 log h 



Figure 1. Upper panel: The SDSS galaxy luminosity function in 
the r-band, K+E corrected to z = 0.1. Data points with error 
bars are taken from Blanton et al. (2003a), as described in §2. The 
shaded region indicates the 95% confidence levels of our Fiducial 
model, as obtained from the corresponding MCMC. Lower panel: 
The corresponding fractional errors, A<I > / < I > . 



P(X\T>) oc exp 



-X 2 (A) 



where X 2 (-M i s a sum °f the following terms 

= XLF + XWp + XESD + Xprior ■ 



(35) 



(36) 



The first three terms quantify the goodness of the fit to 
the data and correspond to the likelihood of the data given 
the parameters, while the last term corresponds to the prior 
information we adopt. The likelihood terms are given by 



3=1 



(37) 



2 

Xwp 



= E t w p( 2fc ) ~ *p( Zfc )l T C wp [wp(%) - Wp(zfc)] ,(38) 



2 

Xesd 



k=l 

6 12 

= EE 

k=l j=l 



AS(fl 3 ,z t )-AS(flj,z t ) 



(39) 



Here X denotes the model prediction for the observable X, 
ax is the corresponding error, w p (zt) is the vector of the 
projected clustering measurements in the fc-th luminosity 
bin, and Cw P is the covariance matrix of these measure- 
ments. 

As detailed in Paper I, our analytical model is accurate 
at the level of 10 (in most cases 5) percent, in reproducing 
the 3-dimensional galaxy-galaxy correlation and the galaxy- 
matter cross-correlation from mock galaxy catalogs. Since 
the differences are not systematic, the accuracy is expected 



to be much better for the projected galaxy-galaxy correla- 
tion function, w p (r p ), and the galaxy-galaxy lensing signal, 
AE(_R). Therefore, we do not account for any systematic 
uncertainty from our modeling in our likelihood estimate. 

For our fiducial model, the set of model parameters, 
A, includes our primary cosmological parameters of inter- 
est, fi m and as, the set of secondary cosmological param- 
eters (fib h 2 , n s , h), the CLF parameters Aclf and a set of 
nuisance parameters (tp,r)). Throughout we adopt uniform, 
non-informative priors on our primary cosmological param- 
eters as well as on all CLF parameters. For the secondary 
cosmological parameters, A scc = (fib h 2 , n B , h), we include 
priors from the WMAP7 analysis. In order to obtain the 
covariance matrix of these parameters (C 8C c), we have used 
the Monte Carlo Markov chains from the WMAP7 analysis 
(kindly provided to us by E. Komatsu) and marginalized 
over all other parameters present in their analysis. Finally, 
for the two nuisance parameters we adopt Gaussian priors, 
as discussed in 53.6. Hence, we have that 



E 



Xprior — [A scc A S cc] C scc [A SC c A SC c] ~t~ 
~ \ 2 

Aj Ai 



0" Xi 



(40) 



where the summation is over the two nuisance parameters 
ip and r). 

We sample the posterior distribution of our model 
parameters given the data using a Monte-Carlo Markov 
chain (MCMC). In particular, we implement the Metropolis- 
Hastings algorithm to construct the MCMC (Metropolis et 
al. 1953; Hastings 1970). At any point in the chain, a trial 
model is generated using a method specified below. The chi- 
squared statistic for the trial model, Xtry> is calculated using 
Equations (36)-(40). This trial model is accepted to be a 
member of the chain with a probability given by 



1-0 if Xtry < Xcur 

CXp[-(Xt 2 ry - X 2 ur)/2] if X?ry > Xcur 



(41) 



where Xcur denotes the x 2 for the current model in the chain. 
We initialize the chain from a random position in our multi- 
dimensional parameter space and obtain a chain of ~ 50, 000 
models. We discard the first 2000 models (the burn-in pe- 
riod) allowing the chain to sample from a more probable 
part of the distribution. We use this chain of models to es- 
timate the confidence levels on the parameters and different 
observables of interest. 

A proper choice of the proposal distribution is very im- 
portant in order to achieve fast convergence and a reasonable 
acceptance rate for the trial models. The posterior distribu- 
tion in a multi-dimensional parameter space, such as the 
one we are dealing with, will have numerous degeneracies 
and in general can be very difficult to sample from. We have 
adopted the following strategy to overcome these difficulties. 
During the first half of the burn-in stage, we chose an in- 
dependent Gaussian proposal distribution for every model 
parameter, as is common for the Metropolis-Hastings al- 
gorithm. Half-way through the burn-in stage, we perform 
a Fisher information matrix analysis at the best fit model 
found thus far. The Fisher information matrix, given by 



d 2 ln£ 

d\id\j 



(42) 
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Figure 2. Left-hand Panel. Projected correlation functions, u> p (r p ), for the six different magnitude bins listed in Table 1, as indicated. 
Solid dots with (1 cr) error bars are the SDSS data taken from Zehavi et al. (2011), while the shaded regions indicate the 95% confidence 
levels obtained from the MCMC of our Fiducial model. We caution that, because of the covariance in the data, which is accounted for 
in the modeling (see §4), the quality of the fit cannot be judged by eye. Right-hand Panel. Excess surface density profiles, AS(H), for 
the six magnitude bins listed in Table 2, as indicated. Solid dots with (1 a) error bars are the SDSS data of Seljak et al. (2005) and 
Mandelbaum et al. (2006). As in the left-hand panels, the shaded regions indicate the 95% confidence levels of our Fiducial model. 



is a N p x N p symmetric matrix, where N p denotes the num- 
ber of parameters in our model, and C oc e~ x / 2 is the 
likelihood. The inverse of the Fisher matrix gives the co- 
variance matrix, C prop , of the posterior constraints on the 
model parameters (see Paper II)**. More importantly, the 
eigenvectors of the covariance matrix are an excellent guide 
to the numerous degeneracies in the posterior distribution, 
and the corresponding eigenvalues set a scale for how wide 
the posterior ought to be in a given direction. Therefore, for 
the second half of the burn-in period, we utilize this informa- 
tion and use a proposal distribution which is a multi-variate 
Gaussian centered at the current value of the parameters 
and with a covariance equal to C prop . In practice, the trial 
model (Atry) can be generated from the current model (A cur ) 
using 

Atry = Acur + C Ax , (43) 



where x is a vector consisting of N p standard normal devi- 
ates, the matrix A is such that AA T = C prop , and c is a 
parameter that we have chosen to achieve an average accep- 
tance rate of ~ 25%. We repeat the Fisher matrix analysis 
once again at the end of the burn-in period (using the best 
fit model found thus far) and use the covariance matrix to 
define our proposal distribution to be used for the MCMC. 
We have found this strategy to be extremely efficient in sam- 
pling our posterior distributions^. 



5 RESULTS 

Having described the data, the model, and the methodology, 
we now turn to our results. In this section we describe the 
cosmological constraints obtained for our Fiducial model, 
whereas §6 discusses the robustness of these results to model 



** The subscript 'prop' indicates that this matrix is used to de- 
scribe the proposal distribution. 



tt A general-purpose python implementation of the MCMC sam- 
pler we have used in our work is available from the authors upon 
request. 
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Figure 3. The 68% and 95% CLs of the joint two-dimensional, 
marginalized posterior distribution for our Fiducial model, ob- 
tained from the simultaneous analysis of the abundance, cluster- 
ing and lensing of galaxies in the SDSS. The green contours show 
the corresponding CLs for the WMAP7 analysis of the CMB (Ko- 
matsu et al. 2011), and are shown for comparison. 



variations. The constraints on galaxy bias, as characterized 
via the CLF, are discussed in §5.3. 

As discussed above, our Fiducial model consists of 16 
free parameters; the two primary cosmological parameters 
of interest, Q m and <78, for which we use uniform, non- 
informative priors, the secondary cosmological parameters 
n s , h and 52b h 2 , for which we use priors from WMAP7 (in- 
cluding their covariance), the 9 CLF parameters that de- 
scribe the halo occupation statistics, also with uniform, non- 
informative priors, and finally the 2 nuisance parameters, i\) 
and r\, for which we adopt Gaussian priors as described in 
§3.6. With a grand total of 182 constraints (32 data points 
for the LF, six bins of 13 data points each for the projected 
correlation function W and six bins of 12 ESD data points), 
this implies 182 — 16 + 5 = 171 degrees of freedom, which is 
the number we have used to compute the reduced \ 2 values 
listed in the final column of Table 3. 

Figs. 1 and 2 compare the predictions of the Fiducial 
model (shaded regions, indicating the 95% confidence levels) 
to the data used to constrain the model (solid dots with error 
bars, indicating the 68 % confidence levels). Fig. 1 shows 
that the model accurately fits the r-band galaxy luminosity 
function. Although most data points agree with the model 
predictions at the 1 a level, the data reveals a few small 
'wiggles' at the faint end that are not reproduced by the 
model, and which contribute dominantly to Xlf> * ne value 
of which is listed in Table 3. 



■ft Although the galaxy-galaxy clustering data points have co- 
variance, we have verified that the covariance matrix for each lu- 
minosity bin has rank equal to 13, and therefore does not reduce 
the number of constraints. 



The left-hand side of Fig. 2 shows the projected correla- 
tion functions, w p (r p ), for six different magnitude bins. We 
caution that, because of the covariance in the data, which 
is accounted for in the modeling (see §4), the quality of the 
fit cannot be judged by eye. However, it is evident from 
the \ 2 values of the best-fit Fiducial model (see Table 3), 
that the total \ 2 is clearly dominated by Xw P - in particular, 
Xw P > 2xesdi even though the projected correlation func- 
tions only have 78/72 ~ 1.08 times as many data points. 
It turns out Xwp is dominated by the contribution from 
the data in the [—20, —21] magnitude bin. Interestingly, this 
bin covers the volume that encloses the Sloan Great Wall 
(SGW), a huge supercluster at z ~ 0.08 and the largest co- 
herent structure detected in the SDSS (Gott et al. 2005). As 
discussed in Zehavi et al. (2011), pruning the data sample 
so as to exclude the SGW region results in a significantly 
reduced clustering strength for galaxies in the [—20,-21] 
magnitude range (i.e., the correlation length is reduced from 
r = 5.46 ±0.15 to 4.82 ±0.23). We return to this issue, and 
its potential impact on our cosmological constraints, in §6.4 
when we discuss the potential impact of sample variance. 

Finally, the right-hand side of Fig. 2 shows the excess 
surface densities, AE(J?), again for six different magnitude 
bins as indicated. The model nicely reproduces the overall 
trends in the data, with only a few data points that fall 
outside the 95% confidence region of the model. Overall, we 
conclude that our Fiducial model is consistent with the 
data at a satisfactory level. In particular, the most impor- 
tant features in the data are nicely reproduced by the model 
and find a natural explanation within the framework of the 
halo model. For example, the fact that brighter galaxies re- 
veal stronger clustering and higher excess surface densities 
is consistent with the common notion that brighter galaxies 
reside in more massive haloes. The lensing signal is directly 
sensitive to this aspect because it probes the matter distri- 
bution around galaxies, whereas the clustering signal is af- 
fected by it only indirectly due to the fact that more massive 
haloes are more strongly clustered than less massive ones 
(e.g., Mo & White 1996). Also, the relatively weak devia- 
tions of w p (r p ) and AE(J?) from pure power-laws typically 
reflect transitions from scales where the signal is dominated 
by different components of the power spectra. Examples are 
the 1-halo to 2-halo transition (e.g., Zehavi et al. 2004) and 
the 1-halo central to 1-halo satellite transition for the excess 
surface densities (e.g., Cacciato et al. 2009). 

5.1 Cosmological Parameters 

Fig. 3 shows the constraints on our two primary cosmolog- 
ical parameters of interest; S2 m and as- The blue contours 
show the 68% and 95% CLs of the joint two-dimensional, 
marginalized posterior distribution obtained from our si- 
multaneous analysis of the abundance, clustering and lens- 
ing of galaxies in the SDSS. The green contours show the 
corresponding CLs for the WMAP7 analysis of the CMB 
(Komatsu et al. 2011), and are shown for comparison. Note 
that our results are in excellent agreement with those from 
WMAP7, strengthening the case for a true concordance cos- 
mology. In particular, our analysis yields S7 m = 0.2781q'q26 
and a 8 = 0.763±%°% (both 95% CL), while the WMAP7 
analysis has Q m = 0.264±g;g|| and cr 8 = 0.801±g;g|| (both 
95% CL). Note also that the degeneracy between Q m and 
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Figure 4. Cosmological constraints for the Fiducial model. Histograms show the marginalized posterior distributions, while the blue 
contour show the 68% and 95% CLs of the joint, two-dimensional marginalized posterior distributions. For comparison, the corresponding 
CLs from WMAP7 are shown as green contours, while the red, solid curves show the marginalized WMAP7 prior distributions used for 
the secondary cosmological parameters, n B , h and f!b h 2 . 



as inherent in our analysis runs perpendicular to that in- 
herent in the CMB data. This indicates that a combined 
analysis will be able to significantly tighten the constraints 
on Q m and as (see also Paper II). Finally, Fig. 3 suggests 
that our constraints are even tighter than those from the 
WMAP7 analysis. However, we emphasize that this is not a 
fair comparison since we have used priors from WMAP7 on 
the secondary cosmological parameters n s , h and f2b h 2 , but 
not on f2 m or as (see Paper II for the case with no priors on 
n s , h and Ob h 2 ) . 

Fig. 4 shows the one-dimensional (histograms) and joint 



two-dimensional (contour plots) marginalized posterior dis- 
tributions on all five cosmological parameters. Solid con- 
tours indicate the 68% and 95% CLs obtained from the 
analysis presented here, while the dotted contours are the 
68% and 95% CLs from the WMAP7 analysis, shown for 
comparison. The strongest parameter degeneracies are be- 
tween fi m and a% (cross-correlation coefficient r — —0.81), 
between fl^h 2 and n B (r = 0.79), and between fi m and h 
(r = —0.74). All other combinations are only weakly corre- 
lated with \r\ < 0.5. 

Overall, there is good agreement between our con- 
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straints and those inferred from the WMAP7 data (see 
also Table 4). However, there is some tension regarding 
the secondary cosmological parameters, which is evident 
from the fact that the posterior and prior distributions 
(indicated by red, solid curves) reveal an offset. This is 
most pronounced for the Hubble parameter h; whereas 
the WMAP7 prior used has h = 0.7101^049 our P oste " 
rior distribution has h = 0.739t^ (both 95% CL). For 
comparison, the revised parallaxes for Cepheid stars by 
van Leeuwen et al. (2007) raises the value for h from the 
HST Key Project from 0.72±0.08kms _1 Mpc -1 (Freedman 
et al. 2001) to 0.76 ± 0.08 km s" 1 Mpc" 1 (68% CL), and 
the Cepheid-based determination of Sandage et al. (2006) 
from 0.62 ± 0.05 kms" 1 Mpc" 1 to 0.70 ± 0.05 kms" 1 Mpc" 1 
(68% CL). Hence, despite some tension with the WMAP7 
based constraints, our posterior distribution for the Hubble 
constant is well within the range of values inferred from 
Cepheids. Interestingly, our constraint on h is correlated 
with the constraints on both fi m and <jg: according to our 
analysis, a relatively low Hubble parameter of h ~ 0.7 im- 
plies f2 m ~ 0.30 and cr 8 ~ 0.73, while for h ~ 0.78 the re- 
sults presented here suggest that $l m ~ 0.25 and erg ~ 0.82. 
Clearly, an improved constraint for the Hubble parameter 
could help to significantly tighten the constraints on both 
f2 m and a$. 

5.2 Nuisance parameters 

As discussed in §3.6, our model contains two nuisance pa- 
rameters: tp, which enters in the description of the scale de- 
pendence of the halo bias, and rj, which sets the normaliza- 
tion of the halo concentration-mass relation. As discussed in 
detail in Paper I, the freedom in r\ also characterizes model 
uncertainties arising from our oversimplifications regarding 
the structure of dark matter haloes. 

The upper panels of Fig. 5 show the posterior distribu- 
tions for rj (left-hand) and tp (right-hand) of our Fiducial 
model. The contours in the other panels show 68% and 95% 
confidence levels of the joint two-dimensional marginalized 
posteriors with fi m and as- Clearly, neither rj nor tp shows 
significant correlation with f2 m or ag (in all cases the cross- 
correlation coefficient \r\ < 0.3). As already discussed in 
Paper II, this is an important result, as it indicates that 
the uncertainties in the scale dependence of the halo bias 
and the oversimplifications regarding the structure of dark 
matter haloes do not have a significant impact on the cos- 
mological constraints presented here. 

The solid, red curves in the upper panels of Fig. 5 re- 
flect the Gaussian priors that we imposed on our model. In 
the case of rj, the posterior distribution (rj = —0.04^']^, 
95% CL) is significantly narrower than the prior distribu- 
tion (rj — O.Otg^, 95% CL), indicating that the prior did 
not have a significant impact on our results. In the case of 
ip, however, the posterior distribution (ip — 1.12*q'i2> 95% 
CL) is clearly offset from the prior (tp = 0.90±g |g, 95% CL) 
to larger values. This might raise concern that a less restric- 
tive prior might have resulted in significantly different cos- 
mological constraints. However, this is not the case for the 
following two reasons. First, as already mentioned above, 
and discussed in more detailed in Paper II, ip is only very 
poorly correlated with the cosmological parameters. Second, 
as discussed in the Appendix, when ip increases the radial 



bias function £(r, z) asymptotes to the empirical fitting func- 
tion of Tinker et al. (2005). In other words, once tp increases 
beyond a certain value, any further increase has zero impact. 
For the best-fit cosmology of our Fiducial model, this crit- 
ical value of ip is ~ 1.5. 



5.3 The Galaxy-Dark Matter Connection 

One of the powerful aspects of the method used here, is 
that the data is used to simultaneusly constrain cosmol- 
ogy and halo occupation statistics. Fig. 6 shows the one- 
dimensional (histograms) and joint two-dimensional (con- 
tour plots) marginalized posterior distributions of model 
Fiducial for the nine CLF parameters that describe the re- 
lation between galaxy luminosity and halo mass. The medi- 
ans and 95% CLs are also listed in Table 5. All nine parame- 
ters are tightly constrained, with tight degeneracies between 
the parameters bo, bi and 62, that describe the normalization 
of the satellite CLF (see Eq. [9]), and between the parame- 
ters Lo, Mi and 71, that describe the relation between halo 
mass and the luminosity of its central galaxy (see Eq. [6]). 

The constraint on the faint-end slope of the satellite 
CLF is a s = -1.18±0.05 (95% CL), in good agreement with 
results obtained from galaxy group catalogues (e.g., Eke et 
al. 2004; Yang et al. 2008). The constraints on the scatter 
in the CLF of central galaxies is a c = 0.157 ± 0.007 (95% 
CL) , which is in excellent agreement with a variety of other 
constraints, from satellite kinematics (More et al. 2009), 
from clustering and lensing (Cacciato et al. 2009; Moster 
et al. 2010), and from galaxy group catalogues (Yang et 
al. 2008). Interestingly, such an amount of scatter is also in 
excellent agreement with predictions from semi-analytical 
models for galaxy formation (e.g., Wang et al. 2006; see also 
More et al. 2009). 

In the left-hand panel of Fig. 7 we compare the con- 
straints on the L c -M relation of our Fiducial model (shaded 
area, indicating the 68% confidence region) to the results ob- 
tained by Yang et al. (2008) from the SDSS galaxy group 
catalogue (symbols with errorbars indicating the 68% con- 
fidence level). Here we have converted the group masses 
listed in their Table 1 to our definition of halo mass and the 
cosmology of the best-fit Fiducial model using the abun- 
dance matching technique described in Eq. (15) of Yang et 
al. (2007). Finally, the right-hand panel of Fig. 7 compares 
the satellite fractions as function of galaxy luminosity ob- 
tained from our Fiducial model (blue shaded area, indicat- 
ing the 68% CL) to constraints obtained by Cooray (2006; 
black, open squares), Yang et al. (2008; red, open circles), 
van den Bosch et al. (2007; green shaded area, indicating 
the 68% CLs), and Tinker et al. (2007; orange shaded area, 
indicating the 68% CL). All these constraints are obtained 
comparing halo occupation models to data from the SDSS 
or 2dFGRS, and are in excellent agreement with each other 
and with the constraints from our Fiducial model. 

Based on all these comparisons, we conclude that the 
constraints on the galaxy-dark matter connection for our 
Fiducial model are in excellent agreement with a wide vari- 
ety of independent constraints. We emphasize that, contrary 
to many other studies, our combined analysis of abundance, 
clustering, and lensing of galaxies also accounts for uncer- 
tainties in cosmological parameters. This strongly supports 
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Table 4. Constraints on Cosmological Parameters 
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The cosmological constraints for the various models discussed in the text (see also Table 3). Values indicate the median and 95% CLs 
of the marginalized posterior distributions obtained from our MCMCs. For comparison, the final row lists the constraints obtained from 
the WMAP7 analysis of the CMB (taken from Komatsu et al. 2011). 
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Figure 5. The shaded histograms in the upper panels show the 
marginalized posterior distributions for the two nuisance param- 
eters r] (left-hand panel), which sets the normalization of the halo 
concentration-mass relation, and if) (right-hand panel), which en- 
ters in the description of the scale dependence of the halo bias. 
The solid, red curves show the corresponding Gaussian prior dis- 
tributions. The contour plots in the middle and lower panels show 
the joint, two-dimensional marginalized posterior distributions of 
Q m and r) (middle left-hand panel), erg and r\ (lower left-hand 
panel), f2 m and ijj (middle right-hand panel) and erg and ip (lower 
right-hand panel). See §5.2 for a detailed discussion. 



that our method of simultaneously constraining cosmology 
and halo occupation statistics yields reliable results. 



6 SYSTEMATICS 

The Fiducial model, discussed in the previous section, re- 
lies on two assumptions regarding satellite galaxies that may 
not be entirely accurate. These concern the radial number 
density distribution of satellites and the Poisson nature of 
the satellite occupation numbers. In this section we gauge 
the impact of these assumptions on our results. In addition, 



we also address the importance of properly accounting for 
the residual redshift space distortions present in the pro- 
jected correlation functions used to constrain the models. 
We do so by running a series of MCMCs in which we in- 
clude small variations with respect to our Fiducial model. 
These models, and their respective x 2 -values for the best-fit 
model in the chain, are listed in Table 3, while the cosmo- 
logical constraints are summarized in Table 4. 

6.1 Radial Number Density Profile of Satellites 

In our Fiducial model we have assumed that the radial 
number density distribution of satellite galaxies in a halo of 
mass M follows an NFW profile with a concentration pa- 
rameter that is identical to that of the dark matter density 
profile, i.e., we assumed that u s (r\M) = ith(rjAZ) (see §3.4). 
This has observational support from a number of studies 
of the radial distribution of satellite galaxies in groups and 
clusters (e.g., Beers & Tonry 1986; Carlberg, Yee & Elling- 
son 1997a; van der Marel et al. 2000; Biviano & Girardi 2003; 
Lin, Mohr & Stanford 2004; van den Bosch et al. 2005b; 
Sheldon et al. 2009). 

However, a number of recent studies have argued that 
the spatial distribution of satellite galaxies is less centrally 
concentrated than the dark matter (e.g., Yang et al. 2005; 
Chen 2008; More et al. 2009). Interestingly, an opposite re- 
sult has been obtained for luminous red galaxies (e.g., Mas- 
jedi et al. 2006; Watson et al. 2010, 2012; Tal, Wake & van 
Dokkum 2012). In fact, there are indications that the ra- 
dial distribution of satellite galaxies appears to have some 
dependence on the color and/or morphologies of the cen- 
tral (e.g., Lorrimer et al. 1994; Sales & Lambas 2005; Lares, 
Lambas & Dominguez 2011; Guo et al. 2012). From a theo- 
retical point of view, one expects the radial distribution of 
satellite galaxies to reflect that of dark matter subhaloes. 
Numerical simulations have shown that subhaloes do not 
follow the same spatial distribution as the dark matter; sub- 
haloes tend to populate preferentially the outskirts of their 
host haloes (e.g., Diemand, Moore & Stadel 2004; Gao et 
al. 2004; Springel et al. 2008). However, the stellar compo- 
nents of satellite galaxies may be more resilient against tidal 
disruption than their dark matter components, such that the 
radial profile of satellite galaxies is not necessarily well rep- 
resented by that of dark matter subhaloes in pure iV-body 
simulations. 
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Figure 6. Constraints on the CLF parameters for the Fiducial model. Histograms show the marginalized posterior distributions, while 
the blue contour show the 68% and 95% CLs of the joint, two-dimaensional marginalized posterior distributions. 



In conclusion, the simple assumption that u B (r\M) = 
Uh(r\M) is likely oversimplified. To address the potential 
impact of this assumption on our results, we now consider a 
model in which we adopt 



Us{r]M) X 



1 + 



TZr, 



(44) 



Here r» is the scale radius of the dark matter density pro- 
file for a halo of mass M, and 1Z is a free parameter. For 
TZ = 1 this profile reduces to that of our Fiducial model, 
i.e., u B (r\M) — Uh(r\M). Hence, 7Z is a parameter that con- 
trols how satellites are concentrated with respect to the dark 
matter. In order to gauge how our results depend on 1Z, and 
whether the data prefers values of 1Z that differ from unity, 



we run a MCMC in which 1Z is a free parameter for which 
we adopt Gaussian prior with mean equal to unity and stan- 
dard deviation, an = 0.25. In what follows we refer to this 
MCMC as our Radial model. 

As is apparent from the x 2 values of the best-fit model 
in the Radial chain (see Table 3), the extra freedom in the 
model does not result in a significantly better fit to the data. 
More importantly, the cosmological constraints are also un- 
affected (see upper left-hand panel of Fig. 8 and Table 4). 
This is also evident from the contour plots in the left-hand 
column of Fig. 9, which show the 68 and 95% confidence in- 
tervals of the joint two-dimensional posterior distributions 
for 7Z and il m (middle row) and for 7Z and erg (bottom row) . 
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Figure 7. Left-hand panel. The luminosity of central galaxies, L c , as function of halo mass, M. The shaded region indicates the 68% 
confidence interval for the Fiducial model, while the open red circles indicate the results of Yang et al. (2008) obtained from the SDSS 
group catalogue of Yang et al. (2007; error bars indicate the 68 % CL), where we properly converted their group masses to our definition of 
halo mass and the cosmology of the best-fit Fiducial model. Right-hand panel. The satellite fractions as a function of r-band magnitude. 
The blue shaded regions indicate the 68 % CL for the Fiducial model. This is compared to the satellite fractions obtained by Cooray 
(2006; black, open squares) from a CLF analysis of the SDSS, by Yang et al. (2008; open red circles) from the SDSS group catalogue 
(error bars indicate the 68 % CL in both cases), by van den Bosch et al. (2007; green shaded area, indicating the 68% CL) and Tinker 
et al. (2007; orange shaded area, indicating the 68% CL), both from halo occupation modeling of galaxies in the 2dFGRS. 



With cross-correlation coefficients \r\ < 0.1, it is clear that 
neither f2 m nor as is significantly influenced by potential 
uncertainty in TZ. Finally, the histogram in the upper left- 
hand panel of Fig. 9 shows the posterior distribution for TZ. 
With TZ = 0.96t° jg (95% CL) it is clear that the posterior 
distribution is narrower than the prior distribution (solid, 
red curve), and that the data prefers values of TZ that are 
close to our fiducial value of unity. We conclude that our re- 
sults are robust to uncertainties in the radial number density 
profile of satellite galaxies, and that the data used here is 
consistent with satellite galaxies being an unbiased tracer of 
the mass distribution within their host haloes (i.e., TZ — 1). 
Interestingly, as mentioned before, this result is supported 
by studies of satellite galaxies in groups and clusters (e.g., 
Beers & Tonry 1986; Carlberg, Yee & Ellingson 1997a; van 
der Marel et al. 2000; Biviano & Girardi 2003; Lin, Mohr 
& Stanford 2004; van den Bosch et al. 2005b; Sheldon et 
al. 2009) but it is somewhat in disagreement with a num- 
ber of studies which indicate that the radial distribution 
of satellite is less centrally concentrated than dark matter 
(e.g., Yang et al. 2005; Chen 2008; More et al. 2009). It is 
worth mentioning that results based on groups might be po- 
tentially affected by the offset of central galaxies (see e.g. 
Skibba et al. 2011). Clearly, more dedicated studies are re- 
quired to further constrain the radial distribution of satellite 
galaxies. 



6.2 Poisson Statistics of Satellite Galaxies 

In our Fiducial model we have assumed that the occupation 
numbers of satellite galaxies obey Poisson statistics, i.e., 



(N S \M) ! 



(45) 



(46) 



P(N a \M) 



such that 

_ (N S (N S -1)\M) 
(N.\M)' 

As shown in Yang et al. (2008), this assumption has strong 
support from galaxy group catalogs. Additional support 
comes from numerical simulations which show that dark 
matter subhaloes (which are believed to host satellite galax- 
ies) also follow Poisson statistics (e.g., Kravtsov et al. 2004). 
However, recently there have been several claims that the 
occupation statistics of subhaloes and/or satellite galaxies 
may deviate slightly, but significantly, from Poisson, i.e., 
Ap ^ 1 (e.g., Porciani, Magliocchetti & Norberg 2004; van 
den Bosch et al. 2005a; Giocoli et al. 2010; Boylan-Kolchin 
et al. 2010; Busha et al. 2011). As shown in Paper I, the 
satellite-satellite term of the 1-halo power spectrum scales 
linearly with _4p, i.e., -Pg^(fe) oc Ap. Hence, any deviation of 
_4p from unity has a direct impact on the projected cor- 
relation functions w p (r p ) on small scales, at least if the 
satellite-satellite term dominates over the central-satellite 
term (which is independent of ,4p). 

In order to quantify how our results depend on uncer- 
tainties related to the exact form of P(N S \M), we run a 
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MCMC in which we include Ap as a free parameter. Moti- 
vated by the empirical results of Yang et al. (2008), which 
suggest that the occupation statistics of satellite galaxies in 
the SDSS are not too different from Poisson, we adopt a 
Gaussian prior for _4p with mean equal to unity and stan- 
dard deviation, oa v = 0.1. In what follows we refer to this 
MCMC as our Poisson model. 

The upper right-hand panel of Fig. 8 shows the 68 and 
95% confidence levels of the joint, two-dimensional marginal- 
ized posterior distribution of £2 m and as for model Poisson 
(yellow contours), compared to that for our Fiducial model 
(blue contours). As is evident, the extra freedom in A p re- 
sults in a best-fit value for as that is slightly lower than for 
the Fiducial model (see Table 4), but the change is only 
marginally significant (i.e., < la). 

The marginalized posterior distribution for Ap, shown 
as the shaded histogram in the upper right-hand panel 
of Fig. 9, has Ap = 0.88loi2 (95% CL), indicating 
that the data prefers a sub-Poisson probability distribu- 
tion P(N S \M). The middle and bottom right-hand panels 
of Fig. 9 show the 68 and 95 % confidence levels of the 
joint, two-dimensional posterior distributions for Ap and f2 m 
(middle row) and for Ap and as (bottom row). Whereas f2 m 
is found to be completely independent of _4 P , there is a weak, 
positive correlation between Ap and as (cross-correlation co- 
efficients r ~ 0.3). Hence, if we would have adopted a less 
restrictive prior on Ap the data most likely would have pre- 
ferred an even smaller value for as- We emphasize, though, 
that our prior has both empirical and theoretical motivation. 
Nevertheless, it is clear that tighter constraints on as require 
a more accurate assessment of the (non)-Poisson nature of 
P(N S \M). 

6.3 Residual Redshift Space Distortions 

In §3.5 we have shown that, when computing the projected 
galaxy correlation functions, we use a modified form of the 
'Kaiser' formalism (Kaiser 1987) to account for the residual 
redshift space distortions (hereafter RRSD) present in the 
data. These RRSD arise as a consequence of the fact that 
the line-of-sight projection was only performed out to a dis- 
tance of r max = 40/i _1 Mpc or 60ft~ 1 Mpc (depending on the 
magnitude bin, see Table 2). In Paper I, we have shown that 
not taking this correction into account results in systematics 
errors that can easily exceed 20 percent on scales above ~ 10 
Mpc. In this section, we show how not correcting for RRSD 
can bias the inference of cosmological parameters such as f2 m 
and as- To this aim, we run a MCMC in which we compute 
the projected galaxy correlation function without correcting 
for RRSD, i.e., we use 

r° ut r dr 

»pM = 2/ g gg (r) Zg— , (47) 

Jr p V 7 " " r P 

with Tout = y r p + r max- As discussed in Paper I, this is 
the standard method used by numerous authors in the past 
(e.g., Zehavi et al. 2004, 2005, 2011; Abazajian et al. 2005; 
Tinker et al. 2005; Zheng et al. 2007, 2009; Yoo et al. 2009). 
In what follows we refer to this MCMC as our noRSSD model. 

The lower left-hand panel of Fig. 8 shows the 68% and 
95% confidence levels of the joint, two-dimensional marginal- 
ized posterior distribution of S7 m and as for model noRSSD 



(yellow contours), compared to that for our Fiducial model 
(blue contours) As is evident, not correcting for the RRSD 
biases as high by ~ 6% and fi m low by about ~ 5% (cf. 
Table 4). Although this systematic bias is only at the ~ 1 a 
level for the current data, its significance will increase when 
using more accurate data and /or data that extends to larger 
radii. Hence, we conclude that, in this era of precision cos- 
mology, properly correcting for RRSD is important. 



6.4 Sample Variance 

Even with a survey the size of SDSS, it is important to ad- 
dress the issue of sample variance. This requires full covari- 
ance matrices for the entire data set (including the covari- 
ance between the LF, WP and ESD data), and requires large 
numbers of realistic mock catalogs. This in turn requires a 
large suite of very large and high resolution iV-body simula- 
tions, which we, unfortunately, do not have available at this 
stage. An alternative is to use jackknife or bootstrap esti- 
mates for the covariance matrix, which can be constructed 
from the data itself. Although these fail to account for sam- 
ple variance on the scale of the survey itself, parameter un- 
certainties estimated from jackknife errors are likely to be 
conservative (Norberg et al. 2011). Unfortunately, since we 
do not have access to the galaxy-galaxy lensing data prod- 
ucts, we cannot construct such a jackknife covariance ma- 
trix. We have recently started a project to perform our own 
analysis of the galaxy-galaxy lensing signal in the SDSS, fol- 
lowing a method similar to that of Mandelbaum et al. (2005, 
2006). Once completed, we will use those data together with 
our own measurements of the clustering and luminosity (or 
stellar mass) functions (see Yang et al. 2012), to construct 
a full covariance matrix for our data set using jackknife or 
bootstrap estimates. 

In this study, however, we have ignored sample variance 
in our analysis, and we caution that our error estimates for 
the model parameters are therefore likely to be underesti- 
mated. However, we can still perform a simple test to gauge 
the potential impact of sample variance on our results. As al- 
ready alluded to in §5, the most important manifestation of 
sample variance in the SDSS appears to be the Sloan Great 
Wall (SGW). In order to get a rough estimate of how the 
presence or absence of the SGW in our data sample impacts 
on the cosmological constraints, we perform the following 
exercise. We construct a MCMC using the same ingredients 
and priors as for our fiducial model, but we exclude the 
w p (r p ) data for the [—20, —21] magnitude bin. As shown by 
Zehavi et al. (2011), this is the magnitude bin that is most 
severely affected by the Sloan Great Wall (SGW). We refer 
to the resulting MCMC as model noSGW. The lower right- 
hand panel of Fig. 8 shows that excluding this part of the 
data lowers the value of as by ~ 5% to as = 0.727^ 049 
(see Table 4), whereas f2 m remains largely unaffected. Note 
that we only removed that subset of the io p (r p ) data that 
is most severely affected by the SGW, but that we did not 
remove the corresponding data from the ESD or LF. The 
cosmological constraints obtained for model noSGW therefore 
have to be interpreted with care. The main purpose of this 
exercise is to demonstrate that sample variance, even in a 
survey the size of SDSS, is still likely to be the dominant 
source of error in the analysis presented here. 
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Figure 8. A comparison of the 68% and 95% CLs of the joint, two-dimensional marginalized posterior distributions of f! m and trg for our 
fiducial model in blue contours and, in yellow contours, model Radial (upper left-hand panel), model Poisson (upper right-hand panel), 
model noRSSD (lower left-hand panel), and model noSGW (lower right-hand panel). In all panels the green contours correspond to the 
constraints from WMAP7, and are shown for comparison. See text for detailed discussion and Table 4 for a listing of the corresponding 
constraints on the cosmological parameters. 



7 CONCLUSIONS 

Galaxies can be observed over cosmological scales, making 
them useful tracers of the underlying dark matter density 
field. Hence, their spatial distribution is a potentially pow- 
erful probe for cosmology. Unfortunately, galaxies are biased 
tracers of the mass distribution, and this so-called galaxy 
bias depends on galaxy properties, on scale, and is even know 



to be stochastic. Because of these complications, galaxy bias 
is generally considered a nuisance when using galaxies to 
constrain cosmology. However, since galaxy bias also con- 
tains a wealth of information regarding galaxy formation, it 
would be highly advantageous if one could simultaneously 
solve for cosmology and galaxy bias. 

Several studies have used galaxy power spectra to simul- 
taneously constrain cosmology and galaxy bias (e.g., Cole et 



© 2008 RAS, MNRAS 000, 1-21 



18 Cacciato et al. 



Table 5. Constraints on CLF Parameters 
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Fiducial 11.24+°^ 9.95l°;°* 3.18±° : » 0.245™ 0.157l°;°°? -1.18l°;° 5 5 -1.17±g;^ 1.53l°;« -0.217±g;°»| 

Radial 11.241°;" 9.95t°;«* 3.17l»;| 0.245^.^ 0.1571°;°°? -1.171°;°° -l.c!4t°;3° USOlg; 8 * -0-212l°;°» 

Poisson 11.23±°;£ 9.95±g : °| 3.30l^ 0.2471°;°" 0.1581°;°°* -^±°o°ol ~^ + -° fs ^t° il -0.2101°-°?? 

noRSSD n.23l«:iJ 9.961°;^ 3.00±°;™ 0.239±°;™ 0.1581°;°°? -1.171°,;°* -1.281 ;g 1.681°;" -0.2621°;?°° 

noSGW 11.021°;!° 9.86l°;° 6 4 4.771°;?° 0.273±g;gl« 0.142l°;°J° -l.0ll«£ -0.64l»;g 1.17l°;« -0.208^;^ 



Constraints on the CLF parameters for the various models discussed in the text (see also Table 3). Values indicate the median and 95% 
CLs of the marginalized posterior distributions obtained from our MCMCs. 
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Figure 9. The shaded histograms in the upper panels show the 
marginalized posterior distributions for the parameter TZ (left- 
hand panel), which controls the radial distribution of satellite 
galaxies (see Eq. [44]), and Ap (right-hand panel), which de- 
scribes deviations of the satellite occupation distribution from a 
Poisson distribution (see Eq. [46]). The solid, red curves show the 
corresponding Gaussian prior distributions used in models Radial 
and Poisson, respectively. The contour plots in the middle and 
lower panels show the joint, two-dimensional marginalized poste- 
rior distributions of Q m and TZ (middle left-hand panel), erg and 
1Z (lower left-hand panel), f2 m and _4p (middle right-hand panel) 
and erg and Ap (lower right-hand panel). See §6.1 and §6.2 for 
details. 



al. 2005; Tegmark et al. 2006; Padmanabhan et al. 2007). 
These studies, however, often resulted in tension regarding 
the cosmological matter density, fi m . This is likely to be a 
consequence of the fact that these studies used a poorly mo- 
tivated, inaccurate fitting function for the scale-dependence 
of the galaxy bias due to non-linearity, causing biased esti- 
mates of the cosmological parameters (e.g., Sanchez & Cole 
2008; Yoo et al. 2009; Hamann et al. 2008). This indicates 
that it is prudent to use a realistic, physically motivated 
model for the galaxy bias. 

In this paper we have simultaneously constrained cos- 
mology and galaxy bias using measurements of galaxy abun- 



dances, galaxy clustering and galaxy-galaxy lensing, all 
taken from the SDSS. We used a realistic, accurate and 
intuitive model for galaxy bias, based on the conditional 
luminosity function (which describes the halo occupation 
statistics as function of galaxy luminosity) combined with 
the halo model (which describes the non-linear matter field 
in terms of its halo building blocks). We improved on pre- 
vious, similar studies (e.g., Seljak et al. 2005; Cacciato et 
al. 2009) by (i) using a more accurate analytical model, (ii) 
using the latest clustering data, (iii) correcting the latter for 
residual redshift space distortions, and (iv) marginalizing 
over uncertainties related to the detailed structure of dark 
matter haloes. 

As detailed in Paper I (van den Bosch et al. 2012), 
our model, which includes scale-dependence of the halo 
bias and a proper treatment of halo exclusion, reproduces 
the 3-dimensional galaxy-galaxy correlation and the galaxy- 
matter cross-correlation (which can be projected to predict 
the observables) with an accuracy better than 10 (in most 
cases 5) percent. Since this uncertainty is not systematic, 
the accuracy is expected to be much better for the pro- 
jected galaxy-galaxy correlation function and the galaxy- 
galaxy lensing signal, which are the observed quantities. Fur- 
thermore, although we believe the halo mass function and 
halo bias function, n(M,z) and bh(M, z), obtained by Tin- 
ker et al. (2008, 2010) to be the most accurate to date, it 
is important to realize that they still can carry uncertain- 
ties that can potentially impact the results on cosmological 
parameters. It is unclear if such uncertainties affect just the 
mass function normalization and not its shape. We intend to 
carry out a proper investigation of this issue in future work. 
Throughout this paper, we restrict ourselves to the n(M, z) 
and bh(M,z) specified above. 

In order to constrain our model parameters, we have 
adopted a full Bayesian framework in which we use Monte 
Carlo Markov Chains (MCMC) to probe the posterior dis- 
tributions. Under the assumption of a spatially flat, vanilla 
ACDM cosmology, and using WMAP7 priors on the sec- 
ondary cosmological parameters, (n s , h, fib h 2 ), we obtain 
that fi m = 0.278t°:oll and a 8 = 0.763l°,;°,t| (both 95% 
CL). These are among the tightest constraints on these two 
cosmological parameters ever obtained, and are perfectly 
consistent with the constraints from WMAP7 (Komatsu et 
al. 2011). We emphasize that no priors were used for either 
fi m or as. Our analysis, which uses a realistic model for 
galaxy bias, therefore demonstrates that the observed dis- 
tribution of galaxies, down to the smallest non-linear scales, 
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is perfectly consistent with the vanilla ACDM cosmology, 
and that the use of a realistic and accurate model for galaxy 
bias removes the tension between the spatial distribution of 
galaxies and the CMB anisotropies that hampered previous 
studies based on galaxy power spectra. 

In order to assess the robustness of our results, we have 
performed a number of tests in which we relaxed some of our 
model assumptions. In particular, we have shown that the 
cosmological constraints are extremely robust to assump- 
tions regarding the radial number density profile of satellite 
galaxies, and that the data is perfectly consistent with satel- 
lite galaxies being an unbiased tracer of the mass distribu- 
tion within their host haloes. Our results are also robust to 
assumptions regarding the second moment of the satellite 
occupation distribution, P(N S \M). In our fiducial model we 
have made the customary assumption that P(N S \M) follows 
a Poisson distribution. Relaxing this assumption lowers the 
best-fit value for a 8 by ~ 2.5% (~ 1 a), to a 8 = 0.744±°;°^, 
while the constraints on fi m remain almost identical. Inter- 
estingly, the data seems to have a slight (marginally signifi- 
cant) preference for a P(N S \M) that is sub-Poisson (i.e., for 
which (N S (N S - 1)|M) < (N S \M) 2 ), consistent with expec- 
tations based on the occupation statistics of (massive) dark 
matter subhaloes (e.g., Kravtsov et al. 2004; Boylan-Kolchin 
et al. 2010). 

We have also shown that not correcting the projected 
correlation functions of Zehavi et al. (2011) for residual red- 
shift space distortions biases as high and f2 m low, both by 
~ 5%. Hence, in this era of precision cosmology, correct- 
ing for these RRSD cannot be neglected. Fortunately, as 
demonstrated in Paper I, RRSD can be corrected for in a 
straightforward and accurate manner using a slightly mod- 
ified version of the Kaiser (1987) formalism. 

The Bayesian analysis presented here has not accounted 
for sample variance. We have performed a simple, naive test 
to gauge the potential impact of sample variance, by exclud- 
ing part of the clustering data that samples the SGW, the 
largest coherent structure detected in the SDSS. This lowers 
the value of as by ~ 5% (~ la) compared to the fiducial 
model, to as = 0. 7271^49, whereas O m remains largely un- 
affected. Hence, we conclude that sample variance may well 
be the dominant source of uncertainty for our results. We 
have recently started to develop our own pipeline to mea- 
sure the galaxy-galaxy lensing signal in the SDSS, which we 
will use to estimate the full covariance for our data set using 
jackknife or bootstrap sampling. This will allow us to better 
address the issue of sample variance for the results presented 
here. 



8 ACKNOWLEDGMENTS 

The work presented in this paper has greatly benefited from 
discussions with Alexie Leauthaud, Rachel Mandelbaum, 
Jeremy Tinker, Risa Wechsler, Idit Zehavi and Zheng Zheng. 
MC acknowledges support from the German-Israeli Founda- 
tion (GIF) 1-895-207.7/2005, and from the Minerva Fellow- 
ship (by Max Planck Gesellschaft) . FvdB acknowledges sup- 
port from the Lady Davis Foundation for a Visiting Profes- 
sorship at Hebrew University. This research was supported 
in part by the National Science Foundation under Grant No. 
NSF PHY11-25915 and NSF PHY-0551142. 



REFERENCES 

Abazajian K., et al., 2005, ApJ, 625, 613 
Abazajian K., et al., 2009, ApJS, 182, 543 
Adelman-McCarthy J.K., et al., 2006, ApJS, 162, 38 
Astier P., et al., 2006, A&A, 447, 31 
Balbi A., et al., 2000, ApJ, 545, 1 

Baldauf T., Smith R.E., Seljak U., Mandelbaum R., 2010, Phys. 

Rev. D., 81, 063531 
Beers T.C., Tonry J.L., 1986, ApJ, 300, 557 
Benjamin J., et al., 2007, MNRAS, 381, 702 
Benson B.A., ct al., 2011, preprint (arXiv:1112.5435) 
Berlind A.A., Weinberg D.H., 2002, ApJ, 575, 587 
Biviano A., Girardi M., 2003, ApJ, 585, 205 
Blanton M.R., et al. 2003a, ApJ, 592, 819 
Blanton M.R., et al. 2003b, AJ, 125, 2348 

Boylan-Kolchin M., Springel V., White S.D.M., Jenkins A., 2010, 

MNRAS, 406, 896 
Bullock J.S., Kolatt T.S., Sigad Y., Somerville R.S., Kravtsov 

A.V., Klypin A.A., Primack J.R., Dekel A., 2001, MNRAS, 

321, 559 

Buries S., Nollett K.M., Turner M.S., 2001, ApJ, 552, 1 

Busha M.T., Wechsler R.H., Behroozi P.S., Gerke B.F., Klypin 

A.A., Primack J.R., 2011, ApJ, 743, 11 
Cacciato M., van den Bosch F.C., More S., Li R., Mo H.J., Yang 

X., 2009, MNRAS, 394, 929 
Cacciato M., Lahav O., van den Bosch F.C., Hockstra H., Dckcl 

A., 2012, preprint (arXiv:1203.2616) 
Carlberg R.G., Yee H.K.C., Ellingson E., 1997, ApJ, 478, 462 
Chen J., 2008, A&A, 384, 347 
Cole S., et al., 2005, MNRAS, 362, 505 
Colless M., et al., 2003, preprint (arXiv:astro-ph/0306581) 
Cooray A., 2006, MNRAS, 365, 842 
Cooray A., Sheth R.K., 2002, Phys. Rep., 372, 1 
Dekel A., Lahav O., 1999, ApJ, 520, 24 
Desjacques V., Nusser A., 2005, MNRAS, 361, 1257 
Diemand J., Moore B., Stadcl J., 2004, MNRAS, 352, 535 
Dunkley J., et al., 2009, ApJS, 180, 306 
Eisenstcin D.J., Hu W., 1998, ApJ, 496, 605 
Eke V.R., Navarro J.F., Steinmetz M., 2001, ApJ, 554, 114 
Eke V.R., ct al. (The 2dFGRS team), 2004, MNRAS, 355, 769 
Freedman W.L., ct al., 2001, ApJ, 553, 47 
Fu L., ct al., 2008, A&A, 479, 9 

Gao L., White S.D.M., Jenkins A., Stoehr F., Springel V., 2004, 

MNRAS, 355, 819 
Giannantonio T., Scranton R., Crittenden R.G., Nichol R.C., 

Boughn S.P., Myers A.D., Richards G.T., 2008, Phys. Rev. 

D., 77, 123520 

Giocoli C, Tormen G., Sheth R.K., van den Bosch F.C., 2010, 

MNRAS, 404, 502 
Gott J.R., Juric M., Schlcgcl D.J., Hoylc F., Vogeley M.S., 

Tegmark M., Bahcall N.A., Brinkmann J., 2005, ApJ, 624, 

463 

Guzik J., Seljak U., 2002, MNRAS, 325, 311 

Guo Q., Cole S., Eke V., Frenk C, 2012, preprint 

(arXiv:1201.1296) 
Guzzo L., ct al., 2000, A&A, 355, 1 

Hamann J., Hannestad S., Melchiorri A., Wong Y.Y., 2008, JCAP, 
07, 017 

Hastings W.K., 1970, Biometrika, 57, 97 

Ho S., Hirata CM., Padmanabhan N., Seljak U., Bahcall N., 2008, 

Phys. Rev. D., 78, 043519 
Jing Y.P., Mo H.J., Borner G., 1998, ApJ, 503, 37 
Kaiser N., 1987, MNRAS, 227, 1 
Komatsu E., et al., 2009, ApJS, 180, 330 
Komatsu E., ct al., 2011, ApJS, 192, 18 

Koopmans L.V.E., Treu T., Fassnacht CD., Blandford R.D., 
Surpi C, 2003, ApJ, 599, 70 



© 2008 RAS, MNRAS 000, 1-21 



20 Cacciato et al. 



Kowalski M., et al., 2008, ApJ, 686, 749 

Kravtsov A.V., Bcrlind A. A., Wechsler R.H., Klypin A. A., 
Gottlobcr S., Allgood B., Primack J.R., 2004, ApJ, 609, 35 

Langc A.E., ct al., 2001, Phys. Rev. D., 63, 042001 

Lares M., Lambas D.G., Dom'nguez M.J., 2011, AJ, 142, 13 

Li R., Mo H.J., Zuhui F., Cacciato M., van den Bosch F.C., Yang 
X., More S., 2009, MMRAS, 394, 1016 

Lin Y.-T., Mohr J. J., Stanford S.A., 2004, ApJ, 610, 745 

Lin H., et al., 2011, preprint (arXiv:1111.6622) 

Lorrimer S.J., Frenk C.S., Smith R.M., White S.D.M., Zaritsky 
D., 1994, MNRAS, 269, 696 

Ma C.-P, Fry J.N., 2000, ApJ, 543, 503 

Maccio A. V., Dutton A. A., van den Bosch F.C., Moore B., Potter 

D., Stadel J., 2007, MNRAS, 378, 55 
Mandelbaum R., et al., 2005, MNRAS, 368, 715 
Mandelbaum R., Seljak U., Kauffmann G., Hirata C.M., 

Brinkmann J., 2006, MNRAS, 368, 715 
Masjedi M., et al., 2006, ApJ, 644, 54 

Mather J.C., Fixsen D.J., Shafer R.A., Mosier C, Wilkinson D.T., 

1999, ApJ, 512, 511 
McDonald P., et al., 2005, ApJ, 635, 761 

Metropolis N., Rosenbluth A.W., Rosenbluth M.N., Teller A.H., 
Teller E., 1953, Chcm. Phys., 21, 1087 

Mo H.J., White S.D.M., 1996, MNRAS, 282, 347 

Mo H.J., van den Bosch F.C., White S.D.M., 2010, Galaxy Forma- 
tion and Evolution, Cambridge University Press, Cambridge 

More S., 2011, ApJ, 741, 19 

More S. van den Bosch F.C., Cacciato M., Mo H.J., Yang X., Li 

R., 2009, MNRAS, 392, 801 
More S. van den Bosch F.C., Cacciato M., Skibba R,. Mo H.J., 

Yang X., 2011, MNRAS, 410, 210 
More S. van den Bosch F.C., Cacciato M., More A., Mo H.J., 

Yang X., 2012a, preprint (arXiv: 1204.0786) 
More S. van den Bosch F.C., Cacciato M., More A., Mo H.J., 

Yang X., 2012b, submitted (Paper II) 
Moster B.P., Somerville R.S., Maulbetsch C, van den Bosch F.C., 

Maccio A.V., Naab T., Oser L., 2010, ApJ, 710, 903 
Navarro J.F., Frenk C.S., White S.D.M., 1997, ApJ, 490, 493 
Netterfield C.B., et al., 2002, ApJ, 571, 604 
Ncyman J., Scott E.L., 1952, ApJ, 116, 144 
Norbcrg P., ct al., 2001, MNRAS, 328, 64 
Norbcrg P., ct al., 2002, MNRAS, 332, 827 

Norbcrg P., Baugh CM., Gaztanaga E., Croton D.J., 2009, MN- 
RAS, 396, 19 

Norbcrg P., Gaztanaga E., Baugh CM., Croton D.J., 2011, MN- 
RAS, 418, 2435 
Oguri M., ct al., 2008, AJ, 135, 512 

O'Meara J.M., Buries S., Prochaska J.X., Prochter G.E., Bern- 
stein R.A., Burgess K.M., 2006, ApJ, 649, L61 

Padmanabhan N., White M., Eiscnstcin D.J., 2007, MNRAS, 376, 
1702 

Padmanabhan N., ct al., 2007, MNRAS, 378, 852 
Park C, Vogeley M.S., Geller M.J., Huchra J. P., 1994, ApJ, 431, 
569 

Peacock J.A., Smith R.E., 2000, MNRAS, 318, 1144 
Percival W.J., et al., 2007, ApJ, 657, 645 
Perlmutter S., et al., 1999, ApJ, 517, 565 

Porciani C, Magliocchetti M., Norberg P., 2004, MNRAS, 355, 
1010 

Pryke C, Halverson N.W., Leitch E.M., Kovac J., Carlstrom J.E., 

Holzapfel W.L., Dragovan M., 2002, ApJ, 568, 46 
Reichardt C.L., et al., 2009, ApJ, 694, 1200 
Reid B.A., Spergel D.N., Bode P., 2009, ApJ, 702, 249 
Riess A.G., ct al., 1998, AJ, 116, 1009 
Rozo E., et al., 2010, ApJ, 708, 645 
Ruhl J.E., et al. 2003, ApJ, 599, 786 
Sales L., Lambas D.G., 2005, MNRAS, 356, 1045 
Sanchez A.G„ Cole S., 2008, MNRAS, 385, 830 



Sandage A., Tammann G.A., Saha A., Reindl B., Macchett F.D., 

Panagia N., 2006, ApJ, 653, 843 
Scoccimarro R., Sheth R.K., Hui L., Jain B., 2001, ApJ, 546, 20 
Scott D., Srednicki M., White M., 1994, ApJ, 421, L5 
Seljak U., 2000, MNRAS, 318, 203 
Seljak U., et al., 2005, Phys. Rev. D, 71, 043511 
Sheldon E.S., et al., 2009, ApJ, 703, 2232 

Skibba R. A., van den Bosch F. C, Yang X., More S., Mo H., 

Fontanot F., 2011, MNRAS, 410, 417 
Smith R.E., et al., 2003, MNRAS, 341, 1311 
Spergel D.N., ct al., 2003, ApJS, 148, 175 
Spergel D.N., ct al., 2007, ApJS, 170, 377 
Springel V., et al., 2008, MNRAS, 391, 1685 
Tal T., Wake D.A., van Dokkum P.G., 2012, ApJ, 751, 5 
Tal T., Wake D.A., van Dokkum P.G., van den Bosch F.C., 

Schneider D.P., Brinkman J., Weaver B.A., 2012, ApJ, 746, 

138 

Tegmark M., Bromley B.C., 1999, ApJ, 518, 69 
Tegmark M., et al., 2006, Phys. Rev. D., 74, 123507 
Tinker J.L., Weinberg D.H., Zheng Z., Zchavi I., 2005, ApJ, 631, 
41 

Tinker J.L., Norberg P., Weinberg D.H., Warren M.S., 2007, ApJ, 
659, 877 

Tinker J.L., Kravtsov A.V., Klypin A., Abazajian K., Warren M., 
Yepes C, Gottlober S., Holz D.E., 2008, ApJ, 688, 709 

Tinker J.L., Robertson B.E., Kravtsov A.V., Klypin A., Warren 
M., Yepes G., Gottlobcr S., 2010, ApJ, 724, 878 

Tinker J.L., ct al., 2012, ApJ, 745, 16 

van den Bosch F.C., Yang X., Mo H.J., 2003a, MNRAS , 340, 771 
van den Bosch F.C., Mo H.J., Yang X., 2003b, MNRAS , 345, 923 
van den Bosch F.C., Tormcn C, Giocoli C, 2005a, MNRAS, 359, 
1029 

van den Bosch F.C., Yang X., Mo H.J., Norbcrg P., 2005b, MN- 
RAS, 356, 1233 

van den Bosch F.C, Yang X., Mo H.J., Weinmann S.M., Maccio 
A.V,, More S., Cacciato M., Skibba R., Kang, X., 2007, MN- 
RAS, 376, 841 

van den Bosch F.C, More S., Cacciato M., Mo H.J., Yang X., 

2012, preprint (arXiv:1206.6890), (Paper I) 
van der Marel R.P., Magorrian J., Carlberg R.G., Yee H.K.C, 

Ellingson E., 2000, AJ , 119, 2038 
van Leeuwen F., Feast M.W., Whitelock P.A., Lancy CD., 2007, 

MNRAS, 379, 723 
Viel M., Weller J., Haehnelt M.G., 2004, MNRAS, 355, L23 
Vikhlinin A., ct al., 2009, ApJ, 692, 1060 

Wang L., Li C, Kauffmann C, Dc Lucia G., 2006, MNRAS, 371, 
537 

Wang Y., Yang X., Mo H.J., van den Bosch F.C, 2007, ApJ, 664 
608 

Watson D.F., Berlind A.A., McBride C.K., Masjedi M., 2010, 
ApJ, 709, 115 

Watson D.F., Berlind A.A., McBride C.K., Hogg D.W., Jiang T., 

2012, ApJ, 749, 83 
Yang X., Mo H.J., van den Bosch F.C, 2003, MNRAS, 339, 1057 
Yang X. Mo H.J., Jing Y.P., van den Bosch F.C, Chu Y., 2004, 

MNRAS, 350, 1153 
Yang X., Mo H.J., van den Bosch F.C, Weinmann S.M., Li C, 

Jing Y.P., 2005, MNRAS, 362, 711 
Yang X., Mo H.J., van den Bosch F.C, Pasquali A., Li C, Bardcn 

M., 2007, ApJ, 671, 153 
Yang X., Mo H.J., van den Bosch F.C, 2008, ApJ, 676, 248 
Yang X., Mo H.J., van den Bosch F.C, Zhang Y., Han J., 2012, 

ApJ, 752, 41 

Yoo J., Tinker J.L., Weinberg D.H., Zheng Z., Katz N., Dave R., 

2006, ApJ, 652, 26 
Yoo J., Weinberg D.H., Tinker J.L., Zheng Z., Warren M.S., 2009, 

ApJ, 698, 967 
York D.G., et al., 2000, AJ, 120, 1579 



© 2008 RAS, MNRAS 000, 1-21 



Cosmological Constraints from Clustering and Lensing 21 



Zchavi I., ct al., 2004, ApJ, 608, 16 
Zehavi I., ct al., 2005, ApJ, 630, 1 
Zehavi I., ct al., 2011, ApJ, 736, 59 

Zhao D.H., Jing Y.P., Mo H.J., Borner G., 2009, ApJ, 707, 354 
Zheng Z., Coil A., Zchavi I., 2007, ApJ, 667, 760 
Zheng Z., Zehavi I., Eisenstein D.J., Weinberg D.H., Jing Y.P., 
2009, ApJ, 707, 554 



depends on cosmology, we include an uncertainty of ~ 15 
percent on ip in our cosmological analysis (i.e., we treat ip as 
a free 'nuisance' parameter, using a Gaussian prior centered 
on ip = 0.9 and with oy, = 0.15). Note that if Eq. (A5) has 
no solution (i.e., when ip is sufficiently large), we set ry = 0, 
which corresponds to simply using the fitting function (A3) 
without modification. 



APPENDIX A: THE CORRELATION 
FUNCTION OF DARK MATTER HALOES 

The function Q(k\M\, Mi, z), which enters in the 2-halo 
terms of the galaxy-galaxy and galaxy-matter power spec- 
tra (Eq. [17]), is related to the power spectrum of haloes 
of masses Mi and M2 at redshift z. Taking account of the 
large scale bias of haloes, its radial dependence (in the quasi- 
linear regime), and the fact that dark matter haloes have 
a finite extent and therefore have a minimum separation 
f'min = r m i n (Mi,M2) (i.e., the so called halo exclusion), the 
two-point correlation function between haloes of masses Mi 
and M2 can be written as 

l + Chh(r,2|Mi,M 2 ) = (Al) 
[1 + 6 h (Mi, z) b b (M 2 ,z) C(r, z) £ mm (r, z)] 0(r - r min ) . 

Here bh(M, z) is the linear halo bias, £(r, z) describes 
the scale dependence of the halo bias in the quasi-linear 
regime, and Q(x) is the Heaviside step function, which as- 
sures that £hh(r, z\Mi, M2) = —1 for r < r m i n in order 
to account for halo exclusion. Throughout we set r m i n = 
max[r2oo(Mi), r2oo(M2)], in accordance with the particular 
halo exclusion adopted by Tinker et al. (2008,2010), whose 
halo mass function and halo bias function we use. As de- 
tailed in Paper I, this implies that 

Q(k\M!,M a ,z) = 

4tt f [l + Z hh (r,z\M u M 2 )] r 2 dr , (A2) 

As shown by Tinker et al. (2005), the radial bias func- 
tion can be written in terms of the non-linear matter corre- 
lation function according to 

_ [1 + 1.17U(V)] L49 (A3) 
C ° ( ' >- [l + 0.69? mm (r,*)p.09- {Ai > 

Here the subscript indicates that this empirical fitting 
function was calibrated using A^-body simulations in which 
the haloes were identified using the friends-of-friends algo- 
rithm. We use a slightly modified version of this calibration 
tuned to reproduce the halo-halo and halo-matter correla- 
tion functions for halos identified by the spherical overden- 
sity algorithm (see Paper 1). This modification only affects 
scales smaller than a characteristic scale, 7y,, such that 

«r,z) = { !! r ^ . (A4) 

The characteristic radius, ry,, is defined by 

log [<o(jy., z) £ mm (rv, z)]=tp. (A5) 

Using high resolution numerical simulations for a ACDM 
cosmology, we found in Paper I that tp = 0.9 can accurately 
describe the halo-halo correlation function on small scales. 
Since we cannot rule out that the calibration parameter ip 
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